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IS.  ABSTRACT 


oratory  experiments  have  established  that  axisyrametric  plasma  accelerator 
using  a  solenoidal  magnetic  'field  can  produce  thrust  by  recycling  the  ambient  material 
of  the  vacuum  tank  in  which  they  were  operated,  indicating  that  an  electromagnetic  accel 
erator  on  a  s/dellite  at  low  altitude  should  be  able  to  ionize  and  accelerate  ;air  strea 
ming  by  it  and  produce  thrust  for  drag  make  up.  This  principle  is  the  basis  of  the 
,Spa<2e /Electric  Ramjet,  (SERJ^Advanced.  Electric..Thru8ter)..  Uti lizat4oriA‘df  solar  cell  arr. 
sto  supply  power  gives  ;the  system  a  theoretical  unlimited  total  impulse;  otherwise  limit 
iby  engine  cathode  and  anode  lifetimes.-  These-lifetime3.can_be.very-long,-however .  , 
^Minimum  system  performance  requirements  for  drag  make  up  in  the  100-300  miles  alti¬ 
tude  range  are  computed. Present  solar  cell  technology  places  a  lower  limit  of  about 
120  miles  altitude^for-lsuccessful  solar-powered  drag  makeup  under  conditions  of  minimum 
drag  by  the'cell* arrays  (alignment  parallel  to  line  of  flight). 

“*^Basic  design  parameters  for  system  components  (anode,  cathode,  magnet  coil,  insulat 
ora)  and  power  requirements  are  calculated.^ aAd  analysis  shews  an  engine  of  less  than 
10  pounds  weight  should  make  up  drag  of  a  1  met er  cylindrical  satellite  for 
;.one  y^ar  at  altitudes  of  120  miles  and  higher. 

Comments  are  made  on  integration  of  the  propulsion  system  into  a  satellite, and  on 
possible  mechanical  and  electromagnetic  interference  effects.  Facility  and  instrumen¬ 
tation  requirements  to  conduct  an  experimental  feasibility  and  performance  program  are 
outlined. 
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FOREWORD 


Recent  satellite  advanced  have  allowed  mission  planners 
to  pursue  mission  goals  that  are  more  ambitious  than  ever  thought 
possible.  Consequently,  there  has  been  a  considerable  impace  on 
satellite  propulsion  systems.  ‘Future  systems  will  be  forced  to 
accomplish  missions  involving  large  total  impulses.  To  reason¬ 
able  achieve  this  performance  level,  the  ambient  atmosphere  should 
be  used  whenever  possible  to  supply  the  thruster  propellant.  Pro¬ 
vided  that  adequate  plasma  can  be  generated  from  the  atmosphere, 
electromagnetic  techniques  can  be  employed  for  drag  make-up  and 
other  low  thrust  maneuvers ,  using  power  collected  from  solar  cell' 
arrays.  The  system  herein  described  has  the  potential  of  pro¬ 
ducing  extremely  large  total  impulses,  extending  the  lifetimes  of 
low  orbit  satellites  by  many  times. 

This  technical  report  has  been  reviewed  and  is  approved. 


WALTER  A.  DETJEN,  Chief 
Engine  Development  Branch 


ACKNOWIEDGSMEMTS 

This  propulsion  oohcspt  tools  definitive  form  while  the 
author  was.  engaged  as  a  visiting  scientist  at:  the  Thermo-Mechanics 
Laboratory  at  Aerospace  Research  Laboratories,  VPAFB.  He  would 
like  to  thank  Mr.  Eric  Soehngen  and  Mr.  Ken  Cramer  of  that  labor¬ 
atory  for  many  valuable  discussions  that  helped  to  crystallize  the 
concept. 


Ever  since  the;  .Initial,  idea  was  presented-  tp:  the_Rocket 
Propulsion  Laboratory*  Major  S.  Baty  has  been  extremely  helpful 
in  suggesting  ideas  and  providing  enthusiastic. -support  for  the 
study.  This  help  is  gratefully  acknowledged. 

Discussions  were,  held  with.  Dr.  Dan^cidin  of  TRW  on  problems 
associated  with  integrating  the  propulsion  system-  into  the  space¬ 
craft.  His  valuable  suggestions  are  acknowledged  with  thanks. 


LIST  OF  FIGURES 
Figure 

2.1  Schematic  of  Engine 

3.1  Anode  Configurations 

3.2  Vapor  Pressure  Curves 

3.3  Cathode  Configuration 

.4.1  F(k^)  vs.  k*-  for  Dipole  Magnetic  Field 

4.2  Lines  of  Constant  Flux  for  Dipole  Magnetic  Field 

7.1  Vehicle  Drag  vs.  Altitude 

7.2  Ratio  of  Anode  Area  to  Solar  Cell  Area  vs.  Altitude 

7.3  Solar  Cell  Area  vs.  Altitude 

7.4  Anode  Area  vs.  Altitude 

7.5  Power  vs.  Altitude 

7.6  Total  Drag  vs.  Altitude 

7.7  Arc  Current  vs.  Altitude 

7.8  Minimum  Discharge  Cross-section  vs.  Altitude 

7.9  Number  of  Turns  of  Magnet  vs.  Altitude 

7.10  Magnet  Mass  vs.  Altitude 

7.11  Magnet  .Surface  Area  vs.  Altitude 

7.12  Magnet  \Temperature  vs.  Altitude 

7.13  Vol.  of  Discharge  Necessary  to  Produce  Ions  for  Exhaust 
Beam  vs.  Altitude 

11.1  Schematic  of  Test  Configuration 

11.2  Performance  of  Aerospace  Chamber  (10V)  for  Space 

Propulsion  Testing  (From  Reference  19) 


v 


LIST  OF  SYMBOLS 


i 


i 

j 

i 


t 


! 


*0 


B 

C 


|e| 


E(k2) 


it 

"D 


T 

F(k) 

i 

I 

k 

k2 

K(k2) 

1 

L 


radius  of  spacecraft  body 
vector  potential 
cross-sectional  area 
magnetic  field  strength 
drag  coefficient 
coefficient  of  thrue.fc. 
voltage  coefficient 
charge  on  the  electron 
electric  field  strength 
elliptic  function 
drag  force 
thrust  force 

grouping  of  elliptic  functions  (see  Eq.4.10) 
induced  current 

electric  current  (discharge  and/or  magnet) 
axial  moment  of  inertia  of  satellite 
gas  constant 

parameter  in  elliptic  function,  - ~ ^ 

elliptic  function 

length  of  anode  ring 

length  of  magnet  wire 

satellite  mass 


vi 


m 

ns 

M 

M* 


n 

N 


P 

P 


<s>. 


*el 


ea 


particle,  mass 
mass  flow  rate 
mass  of  magnet 
Mach  number 

number  density  of  particles 
number  of  turns  on  magnet  coil, 
pressure 
power 


solar  radiant  power  per  unit  area  at  earth's  orbit 
cross  section  for  electron-ion  colisions 
cross  section  for  ionizing  collisions 
radius 

non-dimensional  radius 
resistance  of  electric  circuit 
Reynold's  number 


non-dimensional  satellite  velocity  * 

V  m 

time 

number  of  seconds  in  one  year 

temperature 

gas  velocity 

critical  or  Alfven  velocity 


electric  potential  drop 
ionization  potential 


vii 


axial  coordinate 


number  of  charges  on  ions 
non-dimensional  axial  coordinate 
Various  coefficients,  defined  in  situ 

ionization  rate  parameter 

ratio  of  anode  area  to  solar  ceil  area 

angular  displacement 

magnetic  flux 

angular  displacement 

work  function  of  metal 

electrical  conductivity  of  magnet  wire 

Stefan-Boltzmann  constant 

tsmissivity  of  surface 

viscosity  ?f  gas ;  magnetic  moment 

permeability  of  free  space 

gas  density 

density  of  magnet  wire  material 

fraction  of  atoms  ionized 

surface  reflection  coefficient  for  normal 
momentum  transfer 

surface  reflection  coefficient  for  shear 
momentum  transfer 

solar  cell  efficiency 

period  for  reversal  of  magnetic  field 


viii 


collision  time  for  electrons 
collision  time  for  ions 
cyclotron  frequency  for  electrons 
cyclotron  frequency  for  ions 


TABLE  OF  CONTENTS 

Introduction  1 

Engine  Design  and  Operating  Mechanisms  5 

2.1  Configuration  5 

2.2  Modes  of  Operation  6 

2.2.1  The  Minimum  Power  Hypothesis  7 

2.3  Piasaia  Acceleration  Mechanisms  8 

2.4  Ion  Production  Rate  10 

Electrode  Performance  13 

3.1  The  Anode  13 

3.2  The  Cathode  16 

Magnetic  Field  Considerations  20 

4.1  Characteristics  of  a  Dipole  Magnetic  Field  20 

4.2  Engine  Performance  with  Dipole  Magnetic  Field  22 

4.3  Mass  of  the  Magnet  Coil  24 

4.4  Magnet  Cooling  25 

Vehicle  Drag!  and  Solar  Cell  Performance  29 

Torques  and  Deflections  33 

6.1  Engine  Torque  and  Roll  Angle  33 

6.2  Interaction  of  the  Earth's  Magnetic  Field 

and  Engine  Magnetic  Field  34 

6.3  Perturbations  About  Equilibrium  35 

6.4  Deflections  from  Perpendicular  Magnetic 

Field  Position  36 


x 


7.  Procedure  and  Formulae  for  Evaluating  Engine 

Performance •  Results .  38 

7.1  Satellite  Drag  38 

7.2  Ratio  of  X  .  of  Anode:  Area  to  Solar  Cell  Area  38 

A 

7.3  Solar  Cell  Area  Required  when  Cell  Array 

Parallel  to  Line  of  Flight  39 

7.4  Anode  Area  Required  39 

7.5  Total  Power  Required  by  the  Propulsion  System  39 

7 .6  Total  Drag,  on  the  Vehicle  40 

7.7  Total  Current  in  the  Electric  Discharge  40 

7.8  Minimum  Crdss-Sectionai  Area  for  the  Discharge40 

7.9  Number  of  Turns  on  Magnet  41 

7.10  Mass  of  the  Magnet  41 

7.11  Surface  Area  of  Magnet  41 

7.12  Magnet  Equilibrium  Temperature  42 

7.13  Volume  of  the  Discharge  42 

8.  Discussion  of  the  Calculations  43 

i 

9.  Pulsed  Operation  .44 

10.  Integration  of  the  Advanced  Electric  Thruster 

into  the  Spacecraft  45 

11.  Requirements  for  Wind-Tunnel  Test  46 

11.1  Critical  Engine  Parameters  Requiring 

Investigation  46 

11.1.1  Engine  Components  46 

11.1.2  Performance  Parameters  46 

11.1.2.1  Ionizing  Efficiency  46 

11.1.2.2  Anode  Operation  46 

xl 


11.2  Scaling  47 

11.3  Size  of  the  Testing  Vacuum  Tank  48 

11.4  The  Test  Environment  49 

11.5  Test  Facilities  Available  50 

11.5.1  AEDC  Low  Density  Tunnel  M  50 

11.5.2  AEDC  Aerospace  Chamber  ,(!0V)  51 

11.5.3  AEDC  Aerospace  Environmental 

Chamber  (Mark  I)  51 

11.5.4  NASA  Lewis  Electric  Propulsion 

Facility  51 

11.5.5  Other  Facilities  51 

12 .  Proposed  lExpef imental  Program  52 

12.1  Cathode- Studies  52 

12.2  Investigation  of  Engine  Performance  in 

a  Flow  Field  with  a  Low  Level  of  ionization  52 

12 i 3  Investigation  of  Engine  Performance  in 

Simulated  Flow  Field  53 

12.4  Cost  Estimates  (Submitted  as  Supplement)  -• 
REFERENCES  '  54 

APPENDICES  56 


xii 


1.  Introduction 


'Hie  thrust  to  conduct  the  maneuvers  associated  with  satellite 
reorientation,  orbit  modification  and  drag  make-up  has  up  to  the* present 
been  provided  by  any  one  of  a  number  of  means: 

cold  gas  jets 
hydrazine  jets 
ion  engines 

arc  jets  and  plasma  thrusters  carrying  their  own  propellant 
chemical  rocket  propulsion  systems 
colloid  thrusters. 

The  problems  and  disadvantages  of  these  systems  are  many  and  diverse,  but 
they  have  in  common  the  critical  problem  that  the  total  impulse  available 
is  proportional  to  the  mass  of  propellant  carried  on  the  space  craft. 

Mission  life-times  must  hence  be -programmed  from  the  point  of  view  of 
this  limitation,  rather  than  from  the  mission  requirements.  Utilization 
of  ambient  atmosphere  as  propellant  eliminates  this  problem. 

Studies  have  been  conducted ' to  determine  the  feasibility  of  collect¬ 
ing  air  at  high  altitudes  for  use  in  chemical  or  electric  propul  •'>  on  sys¬ 
tems.  The  most  extensive  investigations  were  those  of  S.T.Demetriades  and; 
colleagues  on  the  Air-Scooping  Orbital  Rocket  (A-SCOR)  (Ref.  1-4)..  Sim¬ 
ilar  studies  were  . carried,  out  elsewhere  (Ref  .5-9) .  These  studies  indi-  t 
cated  that  the  system  was,  feasible  provided  that  a  power  source  of 
considerable  size  (several  Megawatts),  ^was  available  to  .assist  in  the 
collection  process  and  to  heat  the.aif  (thermally  or  electrically)  used 
as  propellant  in  the  thruster  system.  The  lack  of  an  adequate  power  source > 
appears  to  be  the  main  obstacle  to  the  implementation  of  this  system. 

The  ideal  propulsion  system,  for satellites  would  , be  one  for  which  * 
both  power  and  propellant  were available  in  situ..  Rapid  improvements  in 
solar  cell  technology  and  in;  the.  technology  of  axisymmetric  plasma: 
accelerators  over  the  past  ten  years  make  it  desirable  at  this1  time  >t6* 
investigate  the  feasibility  of  such  a  truly  infinite  total  impulse  pro¬ 
pulsion  system.  Obviously;,  unless  extremely,1  large  soier  cell  arrays  are 

•  'i  '  .  -  '  ; 


utilized  (on  the  order  of  square  kilometers) ,  the -power  level  available 
will  be  low  (under  10  kilowatts),  so  that  the  cryogenic  collection  and 
storage  of  propellant--as  .suggested  in  the  above-mentioned  studies  -- 
will  not  be  feasible.  An  alternate  means  of  propellant  handling  is 
obtained  through  the  use  of  a  space-electric  ram  jet,  which  device 
is  herein  investigated.  The  SERJ  utilizes  a  solenoidal  magnetic  field 
and  the  ionization  of  the  ambient  air  in  an. electric  discharge,  as  the 
basic  elements  of  the  operation  of  a  plasma  propulsion  engine,  for  which 
power  is  supplied  by  solar  cell  arrays.  This  system  would  have  an  infinite 
total  impulse. 

1 • i  Axisymmetric  Plasma  Accelerators --Background 

During  the  decade  from  1956  to  1966  both  the  USAF  and  NASA  funded 
extensive  programs  to  develop  low-thrast  electric  propulsion  systems. 
Hampered  by  lack  of  a  definitive  theory  of  operation,,  the  axisymmetric 
plasma  accelerators  lagged  behind  the  development  of  thermal  arc  jets, 
ion  engines  and  colloid  thrusters.  Indeed,  phenomena  were  observed  that 
appeared  to  intimately  couple  the  plasma  engine  operation  and  performance 
with  the  test  environment  (Refs.  10-15')  •  A  great  deal  of  effort  was 
expended  in  attempting  to  decouple  the  engine  from  its  environment  by 
continually  reducing  the  ambient  pressure  into  which  the  engine  was 
exhausted-  Due  to  problems  of  this  nature  the  plasma  thruster,  has  been 
viewed  with  distrust  by  potential’ -users,  and  other  systems  have  been 
given  priority  by  the  development  agencies.  However,  the  unique  character¬ 
istics  of  the  magnetoplasmadynaniic  (MPD)  arc  or  axisymmetric  accelerator 
endow  it  with  some  very  desirable  features  for  low  orbit  missions,  in 
that  it  appears  very  probable  that  it  can  utilize  the  ambient  atmosphere 
as  propellant --in  other  words,  instead  of  attempting  to  decouple  the 
engine  from, its  environment ,  use  the  interaction  as  a  propulsion  mechanism. 
Many  tests  in  various  laboratories  in  this  country  have  established  that 
the.  following  phenomena  occur: 

1.  Thrustvis  produced  by  the  expulsion  of  nigh-velocity 
ions  from  the  engine. 

2.  The  ions  are  produced  by  coliisons  between  atoms  and 
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energetic  electrons  throughout  the  volume  of  the 
discharge . 

3.  The  electrical  discharge  extends  a  considerable 
distance  downstream  of  the  electrodes,  and  hence 
transfers  most  of  the  energy  to  the  gas  far  from 
the  engine.  Viscous  interaction  with  the  engine 
components  can  therefore  be  neglected.  Thermal 
conduction  effects  to  all  engine  components  other 
than  the  cathode  can  also  be  neglected. 

4.  The  solenoidal  magnetic  field  that  is  applied  acts 
as  a  magnetic  nozzle  for  any  plasma  that  is  produced 
by  the  discharge.  That  is,  all  forms  of  electron 

and  ion  energy  (exclusive  of  ionization  and  radiation) 
are  converted  into  axial'  and  radial  ion  velocities 
by  expansion  out  of  the  magnetic  field.  This  indicates 
that  all  the  electrical  power  of  the  discharge  can  be 
converted  into  beam  power  except  for  the  anode  and 
cathode  power  losses,  and  the  power  used  in  ionizing 
the  propellant. 

5.  All  of  the  beam  energy  cannot  be  converted  into  axial 
kinetic  energy  (propulsion),  since  the  ions  must  have 
come  tangential  velocity  to  balance  the  torque  produced 
■bathe  accelerator  by  the  discharge  current  in  crossing 
the  applied  magnetic  field.  Also,  there  will  be  some 
radial  velocity  of  the  ions  due  to  the  shape  of  the 
magnetic  .'ozzle. 

6.  A  large  f ration  of  the  thrust  reaction  occurs  on  the 
magnet,  indicating  that  azimuthal  currents  must  be 
flowing  in  the  plasma  column. 

7.  As  the  mass  flow  injected  through  the  engine  is  reduced^ 
the  discharge  envelope  grows  larger  and  encompasses  more 
mass.  The  high  energy  electrons  in  the  discharge  ion¬ 
ize  some  of  this  gar>  -» nd  use  it  as  propellant.  Eventually 
the  discharge  ionizes  enough  material  to  make  up  a 
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minimum-potential  mode  mass  flow  rate  for  the  engine. 

Under  certain  operating  conditions,  this  total  flow-- 

below  called  the  critical  mass  flow,  or  m  can 
'  .  cr 

be  obtained  from  the  ambient  gas. 

Under  a  wide  range  of  operating  conditions ,  some 
fraction  of  the  discharge  current  is  observed  to  be 
confined  to  a  long  reentrant  filament  that  spins  very 
rapidly  (20khz-500khz)  through  the  gas.  The  filament 
can  be  considered  as  an  ionization  front . 
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2.  Engine  -Design  and  Operating  Mechanisms 
2.1  Configuration 

The  space  electric  'ramjet  (Advanced  Electric  Thruster)  consists 
of  the  following  elements)  (see  Fig. 2.1):  a  component  for  producing. an 
axis ytr*ae trie  magnetio  ;fieiti^,  an  anode  ring^,  a  cathode  assembly 
and  a  power  supply  'Che  magnetic  field  may  be  produced;  by  either 

permanent  magnets  or  by  a  solenoidal  coil.  The  current  to  operate  the 
latter  is  obtained  from  an  auxilliary  power  supply  or  by  the  same ‘used 
to  supply,  the  arc  current.  The  magnets  or  adenoidal  coil  may  be 
placed^ around  the  spacecraft,  or  at  one  or  both  ends.  The  anode  ring 
consists:  of  one  or  more  metallic  rings  plac  >d  concentric  to  the  mag¬ 
netic  field.  Trie  rings  can  be  made  from  any  suitable  material  (e.g. 
brass,  copper,  aluminum,  molybdenum,  iron,  etc.)  depending  on  the 
power  and  discharge  level;  The  size  of  the  anode  ring  is  determined 
primarily^  the  magnitude  of  the  discharge  current,  the  following 
relation  ho lding : 


(2.1). 


where 


I  *  discharge -current 

'»  are«_  of  anode  ring  or  rings 

(n  )  «  electron  density  near  anode  face 
6  A 

|ej »  charge  on  the  electron 
k  »  gas  constant 


m 


mass  of  the  electron. 


Tg  *  electron  temperature  near  anode  face 


anode  attachment  area 
anode,  area-  v* 


Toe  cathode  is  designed  to  produce  electrons  on  the  center  line  of  the 
device  at  a  rate  adequate  to;  carry  the  discharge  current  in  the  vicinity 
of  the  cathode  and  "cathode  jet".  The  cathode  may  consist  of  any  one  of 
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the  following  configurations,  or  any  combination  thereof: 

1.  A  thermionic  emitter  heated  from  some  external  source. 

2.  A  thermionic  emitter  heated  by  energy  from  the  dis¬ 

charge.  This  energy  can  come  from  ion  bombardment  or 
thermal  conduction. 

3.  A  field-emission  type  cathode. 

The  cathode  configuration  should  be  confined  to  as  small  a  cross-sectional 
area  on  the  center  line  (i.e.  concentric  with  the  anode  ring)  as  feasible. 

Electrical  power  to  operate  the  thruster  can  come  from  a  number 
of  sources.  The  most  effective  method  of  obtaining  electric  power  is 
through  the  use  of  solar  cells  and  a  storage  battery,  the  former  either 

on  the' .body  of  the  spacecraft  or  ori  arrays  that  are  deployed  out  some 

distance  from  the  satellite.  A  very  effective  method  of  using  this  type 
of  power  is  to  place  the  satellite  in  an  eccentric  orbit.  During  most 
of  the  orbit  energy  is  being  '.ollected  and  stored  .while  the  thruster  is 
not  operating.  Only  near  pe'igee,  when  the  satellite  is  at  altitudes 
below  400  miles,  is  the  engine  'activated  and  the  resultant  thrust  used 
for  drag  make-up;,  orbit  shifts,  attitude' control,  etc.  Power  can  also 
be  obtained  from  any  other  source  aboard  the  spacecraft,  such  as  fuel 
cells,  batteries,  nuclear  reactors,  etc.  A> promising  new  energy  source 
for  this  specific  application  may  be  the  use  of  high  power  lasers'  to 
beam  power  from  stations  at  the  surface  to  the  spacecraft. 

2.2  Modes  of  Operation  , 

The  thruster  can  be  designed  to  operate  either  on  a  steady  state 
basis,  intermittently,  or  pulsed.  The  type  of  thrust  program  employed 
will  depend^ upon  many  factors,  such  as  power  availability,  possible 
communications  interference,  etc.  One  possible  mode  of  operation  that 
would  be  very  efficient  would  be  operation  of  the  thruster  at  altitudes 
where  the  ambient  ionization  level  is  near  ite  peak,  i.e.  in  the  region 
between  200  and  300  km  altitude.  By  properly  designing  the  magnetic  field 
configuration,  the  discharge  could  encompass  a  large  cross-section  in 
space  and  utilize  or.ly  the  existing  ions  and  electrons.  In  this  manner 
the  thrust  is  produced  without  any  accompanying  i5nization  energy  loss. 
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I 


\ 
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Also,  in  this  mode,  it  will  be  possible  to  operate  at  much  higher  elec¬ 
tric  potentials  and  hence  at  much  higher  specific  impulses. 

2.2.1  The  Minimum  Power  Hypothesis 


For  many  plasma  engines  it  has  been  found  that  the  perfor¬ 
mance  data  can  best  be  correlated  and  understood  if  the  assumption 
is  made  that  the  engine  operates  at  a  minimum  arc  potential  (see 
References  16,17,18).  In  order  to  accomplish  this,  the  power 
input  to  the  beam  and  to  the  production  of  charged  particles  is 
equipartitioried,  i .  e. 


2m 


m-1  (vj  +  vf  +  ...  ) 

m.I  I 


When  this  occurs,  the  exhaust  beam  velocity  is  the  Alfven  speed 


v, 


Al 


/  2  lel  (Vt  + 

V 


+  vf  + .  . . .  ) 


(2.2.1) 


The  mass  flow  rate  in  the  exhaust  beam  adjusts  itself  to  equal  a 
"critical  mass  flow  rate 


This  "adjustment"  can  he  accomplished  in  a  number  of  ways  (see 
Reference  16  ).  For. the  mode  of  operation  under  discussion  in  this 
report,  the  discharge  is  assumed  to  spread  out  in  the  volume  behind 
the  engine  until  it  encompasses  enough  volume  to  ionize  gas  at  a 
rate  sufficient  to  produce  the  critical  mass  flow  rate.  Needless  to 
say,  this  type  of  behavior  has  been  observed  on  many  occasions  and 
often  the  size  of  the  vacuum  tank  was  the  factor  that  limited  the 
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growth  of  the  discharge  volume. 


2.3  Plasma  Acceleration  Mechanisms 

A  very  thorough  discussion  of  these  phenomena  is  given  in  Ref¬ 
erences  13-15,  and  rather  than  reproducing  it  here,  the  relevant  parts 
of  one  of  the  reports  (Ref. 13)  are  appended  to  this  report. 

For  the  specific  case  where  a  deliberate-  attempt  is  made  to  use 
only  the  ambient  gas  as  propellant  and  thus  operate  a  Space  Electric 
Ramjet,  the  following  brief  discussion  outlines  the  engine  mechanisms 
involved . 

There  are  a  number  of  modes  of  operation  for  a  plasma  thruster 
of  this  type,  the  particular  mode  depending  to  some  extent  upon  the 
environment  and  the  configuration.  The  discussion  here  presented  out¬ 
lines  the  generally  accepted  views  of  the  phenomena  involved  and  the 
mechanisms  responsible  for  the  production  of  thrust  when  operating  in  a 
very  low  density  environment  - 

Figure  2.1  illustrates  the  mechanisms  of  the  engine.  Part  of  the 
current  may  be  carried  in  a. "filament"  (which  can  be  many  centimeters  in 
diameter)  that  extends  out  from  the  cathode  and  loops,  back  to  the  anode. 
This  loop  encompasses  all  of  the  magnetic  flux  lines,  inside  the  anode 
ring.  The  "loop"  current,  rotates  rapidly  (200-800khz)  and,  because  of  the 
high  electron  concentration  and  temperature,  within  it,  acts  as  an  ion¬ 
ization  front,  ionizing  the  atoms  of  the  ambient  gas  as  it  spins  through 
them.  The  ions  that  are  produced  are  then  accelerated  (electrostatically) 
toward  the  axis  and  join  with  electrons  from  the  cathode  jet  to  form 
the  exhaust  beam  that  escapes  from  the  magnetic  field  and  produces  a 
net  axial  thrust  and  a  torque  on  the  engine..  It  is  also  possible  to 
consider  that  the  thrust  is  produced  by  j  x  B  forces,  in  which  case  the 
axial  thrust  results,  from  the  interaction  of  the  azimuthal  current  with 
the  radial  component  of  the  applied  magnetic  field. 

The  distance  that  the  current  loop  extends  out  into  the  space  behind 
the  engine;  depends  upon  the  ambient  pressure,  the  applied  magnetic  field 
strength,  and  the  total  current*.  In  some  cases  (very  low  pressures)  it 
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can  be  many  meters  long  and  several  meters  in  diameter  at  its  widest. 

It  is  obvious  that  the  current  can  be  broken  up  into  two 
components :  the  spinning  loop  that  ionizes  the  ambient  gas,  and  the 
ion  current  that  produces  the  thrust.  The  efficiency  of  the  thruster 
depends  upon  the  fraction  of  the  total  current  carried  by  the  ions. 
This  is  usually  between  1/2  and  1/4  of  the  total  current. 

The  engine  produces  a  torque  that  is  equal  to 


BA I(rl  ■  rc> 


(2.2) 


where 


*  magnetic  field  strength  at  anode 
r  *  radius  of  anode  ring 

A 

r  *  radius  of  cathode, 
c 

In  order  to  have  no  torque  on  the  average,  the  magnetic  field  must  be 
reversed  periodically.  Naturally,  use  can  be  made  of  the  torque  for 
spinning  or  despinning  the  vehicle. 

There  are  a  number  of  equivalent  ways  to  calculate  the  thrust 
produced  by  the  engine.  That  most  generally  used  is  a  semi-empirical 


law  that  reads 


thrust  F  *  C_B  Ir 

•  A  A 


(2.3) 


where  *  coefficent  of  thrust,  usually  near  1.  Another  method  is 
to  compute  the  electrostatic  thrust 


where 


r-ff.  lt 

(e|  I  '  nij. 

—  *  mass-charge  ratio  of  the  ions 
*  ion  current 

V  *  anode- to-cathode  voltage 


(2.4) 


B.r.v  C 
A  A  cr  v 


(  with  C  ■  voltage  coefficient  (2.5) 
V  2eV 

cr  raj 

Yj  •*  ionization  potential  of 
ambient  gas  ) . 
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Equations  (2.3),  (2.4),  and  (2-5)  can  be  combined  to  determine  the 
ration  of  ion  current  to  total  electric  current: 


i  -  T 


(2.6) 


This  equation  illustrates  the  very  important  design  criterion  that  the 


product  B.r  should  always  be  kept  low  enough  so  that 

A  A 


|e'VA 

2"lvcr 


<  1 


•  2.7) 


since  otherwise  the  thrust  coefficient  will  drop  to  low  values.  This 
inequality  can  be  written  another  way: 


V-  *=  B.r  v  <  4V_ 
A  A  cr  I 


(2.7a) 


In  order  to  operate  with  high  specific  impulse  and  simultaneously  with 
high  potential  drops,  it  could  be  advantageous  to  multiply  ionize  the 
propellant  so  that 

VI  =  VI^  +  VI2*  +  .  (2.8) 


where 


f  1  \  ch 

Vj  *  1  ionization  potential 
*  2n<*  ionization  optential 


2.4  Ion  Production  Rate 

The  ion  production  is  assumed  to  occur  entirely  in  the  volume 
behind  the  engine  where  the  magnetic  field  is  sufficently  strong  to 
contain  the  discharge.  Assuming  that  the  ions  formed  all  spiral  into 
the  cathode  jet,  the  ion  current :can  be  written  as 

*  Z |e| n  Vol .  (2.9) 

where 
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If  we  assume  the  following,  that 


and 

we  find  that 


and 


n 


1/2 

n/2 


|  kT 


5.6  volts, 


v  *  10°  m/sec 
e 


Vol.» 


I 


2  Z  x  7 .6xl06xl  .6x10  19  x  kT 


Z  $6>2  <  ife> 


m 


(2.15) 
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3.  Electrode  Performance 

The  anode  and  cathode  structures  for  the  SERJ  will  be  in¬ 
vestigated  so  that  the  following  parameters  can  be  established: 

1.  Size  and  weight 

1\  Operating  temperature. 

3.  Power  requirements  and  losses 

4.  Mass  loss  rate  and' .lifetimes. 

In  some  cases,  empirical  information  will  be  necessary  to  accurately 
establish  some  of  the  above  quantities. 

Since  acceleration  of  the  plasma  occurs  by  the  discharge  current 
crossing  magnetic  flux  lines,  it  is  important  that  as  much  of  the  flux 
as  feasible  be  concentrated  between  the  anode  and  cathode  structures. 
Examination  of  the  flux  lines  in  Figure  4.2  indicates  that  the  cathode 
structure  can  probably  be  made  to  encompass  leas,  than  5%  of  .the  flux 
where  <|q  *  4t<jioaNI.  To  accomplish  this,  the  cathode  must  be 
positioned,  slightly  in  front  of  the  plane  of  the  magnet  coil,  on  the 
centerline,  and  have  a  diameter  less  than  0.2a.  The  concentrating  of 
the  flux  lines  around  the  coil  places  severe  restrictions  on  the  size 
and  position  of  the.  anode.  To  have  757.  of  the  flux  between  the  two 
electrode  structures,,  the  anode  must  be  confined  to  the  coil  side  of 
the  0.80  flux  line,  as  might  be  accomplished  by  an  anode  built  as  a 
surface  of  revolution  on  the  *  0.80,  line,  with  a  length  adequate 
to  carry  the  current.  However,  a  short  cylinder  of  diameter  2a  would 
probably  be  adequate. 

Both  electrodes  must  be  placed  in  regions  of  strong  magnetic 
field  to  ensure  that  no  purely  radial  current  flows  between  them.  Also, 
to  prevent  surface  currents  from  flowing  over  the  insulator  between  them, 
the  anode  cylinder  or  ring  should  be  separated  axially  by  a  gap  from 
the  magnet  and  insulator. 

3.1  The  Anode 

A  number  of  possible,  anode  configurations  are  available.  TVo 
that  have  been  much  used  in  laboratory  work  are  shown  in  Figures  (3.1a) 
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and  (3; lb).  The  collection  area  needed  by  the  anode  in  this 
application  would  seem  to  indicate  that  the  ring  configuration 
should  be  used. 

The  anode  collects  electrons  from  the  plasma  to  complete  the 
current  circuit.  There  are  essentially  three^  modes  in  which  it  could 
operate: 

1.  Collection  of  the  current  uniformly  over  the  anode 
surface  with-  no  enhanced  ion  production  near  the  surface 
and  no  sheath  formation. 

2.  Concentration  of  the  current  in  a  filament  and  formation 
of  an  attachment  point.  This  point  would  rotate  rapidly 
around  the  anode  cylinder  due  to  the  j  x  B  forces  from 

the  applied  magnetic  field.  The  current  concentration 
can  occur  only  if  mass  is  available  to  be  ionized  in 
the  filament  near  the  anode  surface.  This  matter  could 
be  ambient  material,  or  mass  eroded  from  the  anode  by 
the  high  heat  f lux  -at  the  attachment  point. 

3.  If,  the  electron  flux  to  the  anode  surface  is  such 

.  that  it  would  carry  more  than  the  total  discharge  cur¬ 
rent,  then  a  sheath  will  form  to  repel!  some  of  the 
electrons. 

Because  of  the  low  ambient  density,  the  low  current  and  the  large 
anode  surface  area,  it  is  extremely  unlikely  that  mode  2  will  ever  occur 
in  this  device.  Further,  mode  3  requires  high  plasma  densities  relative 
to  those  available.  Therefore,  the  analysis  will  be  carried  out  on 
the  basis  of  mode  1  operation. 

The  impingement  rate  .of  electrons  on  the  anode  is 


n  v 

_J§U=.  A 

4  AA 


and  they  carry  a  current 


I  A. 

4  A 


(3.1) 
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where 


|e]  *  electron  (charge  ^  ' 

ng  ■  electron  number  density 
vg  *  electron  thermal  velocity 

«  (  Ma-y* 

v  rf  me  1 
*  anoc?e  *rea> 

The  power  loss  to  the  anode  is.  the  sum  of  the  work  function 
energy  and  enthalpy  carried  by  the  electrons  into  the  metal: 


pa-i<X  +  2!w)  <3'2) 

The  anode  temperature  can  be  estimated  by  equating  the  sum  of  the  power 
dissapated  and  solar  impingent  power,  to  .the  power  radiated  by  the  anode 


m4  _  1  >  dP 


(3.3) 


where  it  is  assumed  that  1/2  of  the  anode  is  exposed  directly  to  the 
solar  radiation. 

The  anode  lifetime  can  be  estimated  by  finding  the  evaporation 
rate  of  the  material  at  its  operating  temperature  T^.  Since  the  vapor 
pressure  is  given  in  atmospheres  in  Figure-3.2,  the  mass  loss  rate  per 
year  of  continuous  operation  is 


where 


so  that 


4  p  t  (£  ) 

P*  y  v  p/A  v 

fl  A 


kgm/yr 


5  2 

p  *  atmospheric  pressure  *  10  newtohs/meter 

m  ■'  *'  '  .'i 


thermal  speed  of  atoms 


t  “  number  of  seconds/year  *  3.15  x  10 1 


1  *  anode  length 


»•*  *  10“  <  |  )A  T 
a  a 


(3.4) 
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For  example,  if  the  anode  is  aluminum,  with  an  area  of  1/4  meters , 
and  runs  at  10Q0°K,  then 


and 


*  883  m/sec 


<  |a>A  '  7  X  10-U 


mloss  ,  1.26  x  1013x  7  x. io-11 
-year  .  '  4 'x  883'* . 

*0.250  kgm/yr  *  6.55  lbs/yr. 

3.2  The  Cathode 

The  conventional  cathode  for  this  type  of  accelerator  would 
be  a  conical  tathode  made  of -thoriated  tungsten;  This  type  has  the 
advantages  of  rugged  construction  and  resistance,  to  breakage,  and ’it 
does  not  require  a  separate  theater,,  since  power  to  liberate  the  electrons 


comes  from  the  electric  discharge  through  thermal  conduction  and  ion 
bombardment .  However,,  such  cathodes  have  been:  operated  only;  at  current 
and  pressure  levels  considerably  higher  than  is  anticipated  here,  and 
there  is  reason  to  doubt  that  the  .same  kind,  of ^performance  can  be  ob- 
tained  at  pressures  less  than  10r  atm  and  currents  of  under  10  amperes. 
Also,  since  the  SERJ  device  would  be  operating  in  air,  the  tip  oxidation 
would  probably  be  rapid. 

For  low  currents  (under  10  amperes)  it  will  probably  be  best  to 
design  the  cathode  with  an  independent  heating  device.  Use  of  a 
thorium  or  oxide  coated  surface  would  minimize  the  temperature  and 
heating  required. 

An  "I/'-type  cathode  configuration*  such  as  that  shown  in  Figure 
3;3,  would  probably  be  optimal.  With  good  insulation,  the  conduction 
power  loss  can  be  kept  to  the  level  of  the  radiation  loss  through  the 
orifice,  hence  the  total  power  loss vwill  be 


where 


K 


(3.5) 
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fsT £=■'■£  temperature  in?  ^theecathode 
-A,,;’*-  catH6de:6rif  1  ce..area 


P^  *  ;pdwer  1088, rfco  the  -cathode. 


Because  of -=the  .pf  jition  pf  theahqde:  relative  to  the  cathode  ,aad  since 
_Jthe  potential  drop  ^vill  be  less^than  lOOyolts,  space  charge  effects, 
will  limit  the  current  of  a  , pure  ^thermionic,  cathode  to  a  few  miiiiam- 
peres'-.  It  will  therefore,  be  necessarjr  to  have  :dn  ion  density  in  the 
cathTodeircayity  (see  Figure  J3, 3).  that  is  comparable  to,  the  electron 
density!.;  These  ddhs  wiTl  consist  of  the  material  used  to  coat  the 
cathode  surface.  The  density^wlll!  be; ;low':eriough  so  that  the  mean 
free  path  of  the  ion-electron  collisions,  will  be  much  larger  than  the 
•cavity  diameter,  so  that  equations  for.  free  particle  flow,  rather  than 
diffusion  equations,  willbe  used  to  describe  the  ion  and  electron  flow 
rates  out  of  the  cavity,  if  both  electrons  and  ions  escape  through  the 
orifice  at  thermal  velocity  the  current  can  be  written  as 


since  the  ion  velocity  is  several  orders  of  magnitude  smaller  than 
the  electron  velocity. 

2 

Cathode  operation  at  a  current  density  of  10  -  20  amps/cm 
through  the  orifice  requires  cathode  temperatures  of  about  2000°K,  so 
the:  electron  thermal  velocity  would  be 


(8  x  1.38'  x  10"23  x  2000l* 
Ve  3il4  x'o:91,x  10-^  J 

*  2v78  x  105  meters/sec. 


The  electron  density  in  the  cavity  is  hence,  from  Equation  (3.6) 
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1.80  x  101S/m3. 


In  order  to  cancel  out  space  charge  effects j  an  ion  .density  comparable 
to  this  value  must  exist  in  the  cavity,  and  this material will-  flow 
out  of  the  cathode  orifice  resulting  in  a  significant  rate  of  loss  of 
material:  - 

V-#  *■  .  0.7, 

For  thorium  and  a  current  of  10  amperes,  the  mass  loss  per  year  would 


^loss  '  1^2 3 2_  x_JJ,6 7  X. x  10  x  3  15  x  in7 

year  1.6  x  10“  19  r  :  x  10  x  x  10 

'*  1.17  kgm/yr  *  2.58  lbs/yr. 

Equation  <3 . 7)  indicates  that  the  mass  loss  rate  is  propor¬ 
tional  to  the  square  rood  of  the  atomic  weight  of  the  cathode  coating 
material.  For  this  reason,  the  feasibility  of  using  a  low  molecular 
weight  substance  rather  than  thorium  should  be  investigated .  For  in¬ 
stance,  if  lithium  were  feasible,  then  the  mass  loss  rate  per  year 
could  be  reduced  to  a  value  of 


/  mloss  v 
^  year  'Li 


1,17  T232 


0.203  kg/yr  -0.448  lb/yr. 


The  momentum  conserve  ion  equation  places  a  further  constraint 
on  the  current  level  at  which  the  L*type  cathode  can  operate.  This 
equation  indicates  that  the  fpl lowing  inequality  mast  be  obeyed: 


(See  Appendix  2.)  From  Equation  (3.6)  an  expression  for  pA  can  be 
obtained  as  follows: 


(3.8) 
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r 


T  *el  aeVe 
1  pA 

,  ve 


8  kT  pA  since  Pj*| 


n  kT. 
e 


Also,  since  'r  ■  ■~r 
~  '  -e  -nSn 

e: 


I  «J£l.  £A 

~\  'rin  v 
e  -e 


Combining  Equations  (3.8)  and  (3.9): 


(3.9) 


or 


%  ; 

8rt  '  lei  Ve 
8?tf?n 


(MO) 

(3.10a) 


For  a  cathode  temperature  of  2000°K,  this  places  a  maximum  limit 
on  the  current : 

T-^  6.28  x  IQ7  X  0.91  K  in'30  x  2.78  x  105 
'  *  ‘  1.6  x  10“19  - - 

<99.3  amperes. 
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4  c  Magnetic  Field  Considerations 

There  are  several  alternative  ways  of  producing  the  required 
magnetic  field.  The  three  most  feasible  systems  would  be 

1.  Permanent  magnets. 

2.  Solenoidal  coil  with  its  own  power  supply. 

3.  Solenoidal  coil  connected  in  series  with  the  discharge.  • 

Since  a  permanent  magnet  would  have  no  power  consumption  ,  it  might 

at  first  glance  appear  optimal.  However,  torques  act  on  the  vehicle 
that  are  proportional  to  the  strength  of  the  magnetic  .field  and  the- 
calculations  in  Section  6  indicate4  that  these  torques  can  seriously 
affect  the  orientation  of  the  vehicle  in  a  relatively  short  time.  The 
only  feasible  method  of  cancelling  these  torques  is  to  operate  the  engine 
for.  equal  times  with  the  two  polarities  of  the  magnetic,  field  direction. 

In  that  it  would  not  appear  practicable  -to  reverse  the  field  of  a  perman¬ 
ent  magnet  configuration  on- a' space  vehicle,  and  because  of  the  weight 
disadvantage  of  permanent  magnets,  such  magnets  would  appear  unsuited 
to  this  application. 

Of  the  two  solenoidal  magnets,  the  seriesrconnected  appears  to 
offer  the  most  advantage.  It  gives  the  electric  discharge  a  positive 
characteristic  so  that  it  can  be  connected  directly  in  series  with  the 
solar  cells,  eliminating  any  power  conditioning  need,  and  It  facilitates 
the  cancellation  of  the  -torques  by  ensuring  equality  of  the  opposite 
torques  (see  Section  6)  in  that  there  is  only  one  current. 

4 . 1  Characteristics  of  a  Dipole  Magnetic  Field; 

The  magnetic  coil  configuration  that  is  most  convenient  and 
probably  near  optimum  in  weight  is  a  short  solenoid,  where  the  coil 
length  is  equal  to  its  thickness.  As  long  as  one  does  not  approach 
too  close  to  the  coil,  it  may  be  considered  a‘  single  turn  of  wire  with 
a  current  of  NI  amperes,  where  N  is  the  total  number  of  turns.  The 
vector  potential  of  the  coil  is  given  by  (c.f.  Jackson,  Classical 
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Elcctrodvruaics .  p  142) 


*0 


4  jiQNIa 

(a2  +  z2  +  r2  +  2ar)* 


(2-k2)  K(k)  -  2E(k) 

-j  (4.1) 


where 

k2  — - - •  l>k^0  (4.2) 

&  +  %■  +  r*  +  2ar 

K(k)  and  E(k)  are  elliptic  integrals  whose,  values  are  tabulated  (as 
in  -Handbook  of  Chemistry  and  Physics)  .  The  magnetic  flux  4  through 
any  area  tir  is  given  by 


t 

1 


f  *  2rfrA0 

- ^oni  it? 

The  magnetic  field  components  are  given  by 
*z  *  r  3?<rV 


^Aq 

“"5T 

-it 


or 


2rfrB 


z  dr 

afcB.  « 


(4.3) 

(4.4) 

(4.5) 

(4.4a) 

(4.5a) 


When  k»  6  (that  is,  for  positions  near  the  s-axis  or  far  from  the 
coil),  the  magnetic  field  parameters  can  be  expressed  simply  as 


*8  '  To”*2"1 

J. 

(4.6) 

(a2  +  z2  +  r2  +  2ar)3/2 

2(a2  +  z2)  -  r2  +  ar 

(4.7) 

(a2  +  z2  +  r2  +  2ar)5/2 
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3zr 


a 

Br  "  NI  ,  2  ,  2.2.,  .5/2 

(a  +  z  +  r  +  2ar) 

(4.8) 

2  2x11:2 

+  ■  JV4*  1,1  (a2  +  1  +  r2  +  2ar)3/2 

• 

(4.9) 

The  vector  potential  and  flux  can  be  evaluated  for  the  region  around 
the  wire,  numerically.  Defining 


p(k)  „  (2  -  K  )K(k)  -  2E(k) 


(4.10) 


then 


*  4i^ioNI  Var  F(k)  •  (4.3a) 

The  function  F(k)  is  plotted  in  Figure  4.1.  For  values  of  k<0.5, 
the  function  F(k)  can  be  expressed  accurately  by  a  series 

F (k)  V?  rt  (  |2)  3/2  f  1  +  6(  )  +  ~(  y  )2  + 

k2  3  1  (4-11) 

245  (  F)  3  . } 


The  coordinates  of  constant  flux  lines:  are  shown  in  Table  4.1.  Lines 
of  constant  flux  are  plotted-  in  Figure  4.2. 


4 .2  Engine  Performance  with  a  Dipole  Magnetic  Field 


At  this  point,  it  is  desirable,  to  express  the  engine,  performance 
in  terms  of  the  magnetic  dipole  parameters.  Using  the  results  from 
Appendix  2,  the  torque  on  the  engine  can  be  expressed  as 


T  = 


1  ^AC 


(4.12) 


where 


I  »  current  of  discharge 

and  (^  •  »  minimum  flux  enclosed  between  the  anode  and  cathode 

surfaces'; 

Equation  (4;12)  is  a  generalization  of  Equation  (2.2),  A  generalized 
expression  for  the  thrust  to  replace  Equation  (2.3)  can  also  be  written 
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in  terms  of  the  magnetic  flux,  as  follows: 


1  ♦. 


AC 


T  it'  r, 


or  alternatively 


1  i 


AC 


■n  r. 


(4.13) 


(4.13a) 


where  r_  m  reference  radius.  If  a  cylindrical  anode  is  used,  the 

iC 

minimum  value  of  r^  would  be  the  anode  radius  r^.  The  maximum  radius 
iy  could  hove  would  be  the  radius  at  which  the  flux  line  through 
the  rear  of ;  the  anode  is  perpendicular  to  the  coil  centerline.  This 
value  can  be  determined  from  Table  4.1.  In  general,  the  reference 
radius  will  increase  as  the  ambient  density  decreases,  hence  less 
thrust  is  produced  in  a  low  density  environment  than  would  be  produced 
in  one  of  higher  density.  For  this  reason  we  shall  use  the  most 
pessimistic  value  for  (rp)max““^*e*  t*le  that  gives  the  minimum 

thrust  for  a  given  amount  of  magnetic  flux,  between  the  cathode 
and  anode. 

In  order  to  estimate  the  size  of  the  magnet  required,  the 
following  assumptions  are -made: 

1.  The  magnet  and  engine  are  run  in  series  so  that  the 
magnet  current  and  discharge  curre'ut  are  the  same. 

2.  The  anode  is  a  cylinder  of  1  meter  diameter  and  1/4 
meter  long,  and  insulated  on  the  outside. 

3.  The  magnet  coil  has  an  average  diameter  of  1  meter. 

4.  The  front  end  of  the  anode  is  placed  near  the  magnet 
coil. 

5.  The  flux  encompassed  by  the  cathode  is  negligible. 

The  minimum  magnetic  flux  passing  through  the  anode  can  now  be  evaluated, 
using  Equation  (4.3a): 


Iac  - 4  *  N  ra«^A  P<V 


R. 


rA 

a 


-  1 


(4.14) 


where 

and 
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-0.9412 


kA“(1  +  V2+2 * 

2a 

Z.  =  —  =0.50. 

A  a 

Using  the  tables  for  evaluation  -,.F(k^)  is  found  to  be  0.888.  Hence 


*  0.888  x  4i^0NIa 


(4.14a) 


Following  the  flux  line  that  passes  through  the  rear  of  the  anode,  it 
is  seen  to  become  perpendicular  to  the  magnet  centerline  at  approx¬ 
imately. 

2^1  =  0.59 
«  1.38. 

The  thrust  equation  (4.17a)  can  now  be  rewritten  as 


u  I  x  4t1uq  NIa  x  0.888 
=  r  1.38nra 


F  =  2.57  u  NI2 

_ 111 J-° _ 


(4.15) 


F  -  3.23  x  lO'^NI2. 

This  equation  can  be  used  to  evaluate  the  number  of  turns  N  required 
on  the  magnet. 


4.3  Mass  of  the  Magnet 

The  mass  of  the  magnet  can  now  be  estimated.  Using  Ohm's  Law 

L 

where 

=  electrical  conductivity  of  magnet  wire 
Ay  =  cross-sectional  area  of  magnet  wire 
V  =  voltage  drop  across  the  magnet  coil 
L  -  length  of  magnet  wire. 
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The  mass  M  of  the  magnet  coil  can  be  written  as 


M  -  o-  L 


(4.17) 


Eliminating  A  between  Equations  (4.16)  and  (4.17)  gives 

m£u  l2  i2 

°W  PM 


(4.18) 


If  the  magnet  is  wound  with  a  diameter  close  to  that  of  the  vehicle 


L  «-rfDN 


and  the  magnet  mass  can  be  written  as 


fu  ftfflsi)2 

°W  PM 


(4.18a) 


The  parameter  that  determines  the  best  material  to  use  for 
the  magnet  is  the  ratio  of  the  density  to  the  electrical  conductivity. 
Values  of  this  parameter  for  verious  metals  are  presented  in  Table  4.2. 
The  values  indicate  that  the  lightest  magnet  would  be  made  of  sodium, 
which  metal,  if  used,  would  require  enclosure  by  a  tube  (e.g.  of  stain¬ 
less  steel)  to  prevent  breaks  in  the  circuit  due  to  melting  or  reacting 
of  the  sodium. 

4.4  Magnet  Cooling 

The  magnet  will  be  radiation  cooled.  The  temperature  at  which 
the  coil  must  be  run  in  order  to  handle  the  power  can  be  computed  from 
Stefan's  Law: 


where 


P  «  A  e  o'r 
M  s  R 


(4.19) 


power  radiated 


surface  area  of  wire 


o'*  Stefan-Boltzmann  constant 


****** 


f 


T  “  wire  temperature 
eR=  emissivity  of  wire  surface. 
iThe  surface  area  of  the  wire  is  given  by 


•  AT«aii 

o  W 


where  r^  is  the  wire  radius.  But 


so  that 


M  (2rfaI^)4 


Hence  we  have 


r  =  (.  — --  t  JL 

-  '  2rtaN  '  v  ) 

Jv 


w 


(4.20) 


(4.21) 


M  (  £-  )%  (2rt.»)‘> 
J  w 


2rf..( 


2aMN  % 

f w 


(4.21a) 


If  the  coil  is  wound  in  a  single  spir.1,  the  auri.ce  .red  thet  c.n 
radiate  heat  will  be  approximately  h*lf  the , total  wire  surface.  When 
there  .re  many  turns,  the  coll  must  be  ound  in  leyeri,  and  if  it  is 
wound  so  thit  its  thickness  a«d  width  are  equal,  the  surface  are.  will 


'Vein  -  8  <  ~  >* 

J" 


(4.21b) 
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5.  Vehicle  Drag  and  Solar  Cell  Performance 

The  drag  force  that  must  be  balanced  by  the  engine  thrust  con¬ 
sists  of  the  drag  on  the  space  vehicle  itself  and  that  on  the  solar 
cells  and  other  engine  components  which  supply  the  power  to  the  engine. 
The  formulae  needed  to  estimate  the  drag  are  given  below. 

The  pressure  on  a  surface  due  to  molecular  impingement  is 

given  by 


+  VXS'>2 

2,  T  ' 


,  \  pv2  I"  (2  -  <r')S' 
v+|  (2  -  <r  ')  (  S'  2  + 1  )  +  fj*/(  ^  )k  sj 


(5.1) 


i  +  erf (S' 


} 


The  shear  stress  on' the  surface  is 


«  |  .^2r08f  (e"(S,)2  *  **  +  erf(S*>  )J  (5.2) 

where  . 

s2  -  2C 

■  2kT 


S'  -  S  siny 


v  *  satellite  velocity 


“  satellite  surface  temperature 
y  *  angle  between  surface  and  incident  beam 


gas  temperature 


2  f£'  -<S’)2dS' 

erf (S')*  J  *• 


o 

9“  o'1*  *  surface  reflection  coefficient  for  normal  momentum 
transfer,  shear  momentum  transfer 
As  reference  drag  force  for  the  vehicle  we  will  assume  that  the 
vehicle  shape  is  that  of  a  cylinder  one  meter  in  diameter  and  one  meter 
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in  length,  and  that  it  is  aligned  so  that  the  cylinder  axis  is 
parallel  to  direction  of  flight.  Formulae  (5.1)  and  (5.2)  reduce, 
under  these  assumptions , “to 

P  »  2  (2  -  V')  \f  (5.1a) 

”  %Ts2  f  \  <5‘2*> 

The  total  drag  force  on  the  cylinder  can  therefore  be  computed  and 
is  given  by 

Fb  -  {2(2  -  <r  •)£  D2  +  (5.3) 

•  (2(2  -  o' ')  f  O2  5  f  V2  (5,3.) 

This  allows  us  to  compute  the  drag  coefficient  C^: 

Cd»2<2  -<r-'>  +£sf  (5.4) 

■  2<2  -  «•««)  4 
■  2.23  <•  0.25 
-  2.48. 

If  solar  cell  arrays  are  used  to  collect  power  to  run  the 
propulsion  system,  their  drag  must  be  computed.  This  drag  is  strongly 
dependent  on  the  orientation  of  the  array  relative  to  the  direction  of 
satellite  motion.  If  they  are  perpendicular  to  the  motion  they  will 
have  a  drag  coefficients^  «  2.23,  while  alignment  parallel  to  the 
motion  gives  CDsc  »  0.127.  The  type' of  orbit  that  is  used  will  deter¬ 
mine  the  angle  of  the  solar  cells,  so  that  no  attempt  to  find  optimum 
angles  is  made  here.  Rather  performance  capability  for  both  normal  and 
parallel  orientation  will  be  presented. 

The  power  collected  from  the  solar  cells  can  be  written  as 


where 


(  ^  )'  '■  power  density  of-  solar  energy 

i^sc  *  efficiency  of  solar  cells- 

A  *  cross-sectional  area  of  solar  cells., 
sc 

Assuming  equipartition  of  the  power  between  losses- and  beam 
power,  the  power  used  by  the  engine  can  be  written'  as 


PM+22A  +  ?E 
F2 

pm  +  T 


(5.6) 


where-  P£  *  electrode  power 

Pw  ■  magnet  power 

M 

F  *  thrust 
m  ■  mass  flow  rate. 

For  preliminary  estimates  we  can  let  the  sum  of  the  magnet  and 
electrode  power  be  equal  to  the  engine  power,  so  that  the  overall 
efficiency  of  the  engine,  is  25%.  Later,  the  weight  of  the  magnet 
required:  to  permit  this  performance  will  be  estimated. 

The  thrust  F  can  be  written  as 

F  *  m  v  «-m  (v  —  y  )  (5k7) 

cr.  '  e  v 

Equating  this  to  the  total  drag  gives  the  power  balance 

<  4*  >.VA.o  ’  2  V«  x  Dt*8 

'  2  vcr  2  /  \  { <VA  +  <VscAsc  (5 ,8) 
+  <CdVA} 
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mmsssissm. 


The  drag  of  the  anode  has  been  added.  Since  it  is  the  only  other 
major  component  of  any  size  required  for  the;  propulsion  system,  the 
total  drag  of  the  vehicle  and  propulsion  system  has  been  found. 

The- anode  can  be  designed  in  several  ways,,  one  an  annular 
plate  with  the  outside  diameter  approximately  equal  to  that  of  -the 
vehicle  (see  Figure  3.1a).  This  design. has  the  disadvantage  that  the 
anode  surface  area  is  limited  and  may  not  be  large  enough  to  collect 
all  of  the  current..  Ah  alternative  design  is  a  thin  cylinder  with  a 
diameter  approximately  equal  to  that  of  the  vehicle  (  Figure  3.1b) , 
a  configuration  permitting  large  surface,  areas  with  minimum  drag. 

For  such  anode  configuration  the  surface  area  needed  to  collect 
the  current,  can  be  expressed  in  terms  of  the  solar  cell  surface 
area  A  as  follows: 

SC 

I_  *  m  *  ioh;  current  (5.9) 

l  m 

I  aa  '»  total  current  (5.10) 

Combining  Equations  (5.7),  (5.9),  and  (5.10)  gives 

i  _ ft _  1  ZJeJ  JL 

A  ie|neye  Ij  m  ^ ' 

-JUL.  I  JL  . 

"Ve  ZI  vcr  .  (5‘U> 

Using  Equation  (5.8),  this  can  be  written  as 
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Defining 


gives 


2  .  I  dP  . 

?pv  v,  IT  -  dA  '  s  T 
J  e  cr  I 


sc 


4*: 

A  sc 


(5.13) 


(544) 


^  can  be  evaluated  as  a  function  of  altitudeonce  values  are  assumed 
for  e  and  l/r_._  In  the  following  calculations  we  will  assume 


•2 

2. 


With  the  anode  cylinder  aligned  parallel  to  thedirection  of  motion* 
the  drag  coefficient  of  the  anodes  (C^)^  will  be  the  same  as  that  for 
the  solar  cells.  This  permits  calculation  of  the  solar  cell  area 
using  Equation  (5.8): 


Vv 


sc 


1  dA  £:•: 


Jail  s.1 


(5.8a) 


2v 


cr 


If 


(Cd^(1^a> 


When  the  solar  cells  are  normal  to  the  line  off  light,  the  maximum 
density  at  which  the  solar  cells  can  supply  adequate  energy  for  drag 
make  up  is 


/Max 


i£  \  .  “Has 

dA-  »  2  vcr 

2.23  v2  fi 


(5.15) 


The  solar  energy  constant  Is 

(4?  )  «  1300  watts/m2 

dA  s 

and  for  the  other  quantities  we  can  take  as  approximate  values 


V0-10 

vy  ■  8000  meters/sec 
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and  obtain 


v.  *  1.6  x  104  meters/sec 
cr 

,  1.3x  103  x  10*1  1 

jW  2  X  1.6  X  1C*  1.11  X  o.64  x  10$ 

*5.7  x  10  kgm/ra^. 

This -.corresponds  to  an  altitude  of  approximately  300  km  (186  mi).  For 
a  critical  velocity  of  2.4  x  104  meters/sec  (  doubly  ionized  atoms) 
the'  solar  cells  can  supply  power  to  make  up  their  own  drag  only  at 
higher  altitudes  (above  210  mi)  if  the  solar  panels  are  normal  to  line 
of  flight. 


6.  Torques  and-  Deflections 
6 • 1  Engine  .Torque  and  'Roll  Angle 


The  engine's  electromagnetic  torque  will  tend  to  spin  the 
vehicle,  and  if  for  a  particular  satellite  this  spin  is  considered 
undesirable,,  periodic  magnet  current  reversal  will  be  necessary  to  keep 
the  angular  rotation  or  roll  §  about  the  vehicle,  axis  to  within  prescribed 
limits.  The  equation  of  motion  for  this  deflection  is 


d2j 

-  T.  (6.1) 

2  *  dt2  * 

where  J^“  mv  ~~  * 8  the  axial  moment  of  inertia  of  the  vehicle  la  terms 
of  its  outer  radius  a,  and  T^  «  Ihz  from  Equation  (4 . 12) .  Integrating 
equation  (6.1)  and  specifying  a  maximum  permitted  angle  of  roll  for 
the  vehicle,  the  period  'X  for  reversing  the  magnet  current  can  be  estima¬ 
ted 

X  -  (6.2) 


From  equation  (3.14a) 

so  that 


As  an  example, 


let 


(6.2a) 


T  1 

0  ttt  radians  (about  1  degree) 

rtu  ou 

*  100'  kgm 


and 

These  give 


a  *  1/2  meter 
NI  *  100  ampere- turns 
I  »  10  amperes . 


rjs  f  100  x  .  5  x  1 

1  [1.776x1.25x1c'6 


x  60  x  10 


1 


<  1C 

—  15  sec. 
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6.2  Interaction  of  the  Barth's  Magnetic  frield  and  Engine  Magnetic  Field 

There  will  be  two  torques  on  the  vehicle  due  to  the  engine  magnet 
interacting  with  the  earth's  magnetic  field.  When  the  coil  has  an  applied 
current  flowing  through  it  there  will  be  a  torque  tending  to  line  up  the 
two  fields.  There  will  also  be  a  "  damping''  torque  on  the  vehicle  due 
to  the  induced  current  in  the  coil  caused  by  the  change  in  flux  through 
the  c oil  as  the  vehicle  rotates  in  the  earth's  magnetic  field.. 

The  torque  caused  by  the  current  through  the  magnet  can  be  easily 
evaluated:  it  is  the  product  of  the  magnetic  moment  of  the  coil  (^t) 
and  the  terrestrial  field  strength; '(B)  times  the  sine  of  the  angle  be¬ 
tween  them  (  0  ) , 


where 


B  sin  0 
e 


u  -  NI  A 
/  a  c 

N  *  number  of  turns  on  coil 

A  *  area  of  coil, 
c 


(6.3) 

(6.4) 


If  the  fields  are  antiparallel,  the  negative  sign  should  be  used  in  (6.3), 
and  if  parallel,,  the  positive  sign. 

The  damping  force  can  be  evaluated  by  using  the  induction 


equation'  , 

* 

1R*-  -H 

.  (6.5) 

where 

is  the  wire  resistance 

with 

w  w 

i  *  induced  current: 

2 

4  *  flux  through  coil  *  -n'  a  B  cos0 
Aw  *  cross-section  area  of  magnetic  wire 
■  conductivity  of  magnet  wire 
a  *  radius  of  coil. 

The  induced  current  can  therefore  be  expressed  as 

2na  e  dt 


(6.6) 


A  o^aB  siri0 
w  w  e  d0 
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T  «  +  i  rU2m 
2  —  e 

ilk 

"  a  h3e8in0  dF 


(6.7) 


A  differential  equation  for  the  motions  of  the  vehicle  due  to  these: 
torques  can  be  set  up: 


J  +  ~A  o-a3NB2sin0  77  +  NItk2B  sin9  -  0 
dt2  2w  w  e  dt  e 

W2 

_ a  _  -i.  _ c  i _ r  *.i 


(6.8) 


where.  J  »^ii8  is  the  moment  of  inertia  of  the  vehicle  about  one  end, 
assuming  vehicle  of  length  2a  i  Equation  (6.8)  can  ;be  rewritten  as 


V  2 

i. *  +  &.WV  a»  +  i2sr  MB.sl.i9 

'  2  19  mv  81  dt  19  m  ' 

at 


(6.8a) 


which  is  ah  equation  very  similar  to  that  describing  the  motion  of  a 

simple  pendulum  under  gravity.  Exact  solutions  are  not  possible,  but 

approximate  solutions  can  be  obtained.  The  case  of  small  deflections 

from  either  the  stable  position  (0*0) ,  where  the  fields  are  aligned 

'ff 

antiparallel, >  the  orthogonal  position  (0*j),  ind  the  unstable  equilibrium 
position  (0"ri)  are  considered  below.  The.  first  and  third  of  these  cases 
would  correspond  approximately  to  a  satellite  in  a  polar  orbit,  and  the 
second  to  one  in  an  equatorial  orbit . 


6.3  Perturbations  About  Equilibrium 

In-  this  case  we  make  the  approximation  that  sin0  =  0  and  neglect 
the  damping  term.  Then 


iLf  +  A2S  Sfia-e^o 

.,2  19  tv 


(6.8b) 


and  so 


6  ■  B1  *  +  e2  1  (6-9) 


representing  oscillations  with  a  period  of 


m. 


12n  KIBe 


(6.10). 


As  an  example,  let 

mv=  100kg  (220  lb) 

NI  =  100  ampere- turns 
-4 

Bg  =  10  Tesla  (1  gauss) . 

These  values  give 

=  445  seconds . 

6.4  Deflections  from  Perpendicular  Magnetic  Field!  Position 

In.  this  case,  sin9  =  1.  If  the  damping  is  neglected,  then  the 
solution  to  the  Equation  (6.8a)  is  found  tp  be 


6n  NIBp  2 
"  19  iv 


(6.11) 


The  time  2^  required  for  a  deflection  of  an  angle  away  from  the 
orthogonal  position  can  now  be  computed  to  be 


(6.12) 


Using  the  same  values  as  .given  previously  arid  permittirig  a  deflection 

of  9  “  -rj:  (about  1  degree) ,  the  time becomes 
m  ou.  '  o 


r0  Y 


19  x  100 


100  x  60  x  lO-4  x  18.85 


seconds 


=  12.9  seconds. 

For  a  vehicle  in  an  equatorial  orbit,  the  magnetic  field  should  hence  be 
reversed  about  once  every  10-20  sec  in  order  to  keep  this  deflection 
lower  than  1-2  degrees  from  the  desired  orientation. 

6 -5  Deflections  from  the  Unstable  Equilibrium  Position 

Here,  0  and  any  deflection  generates  a  force  tending  to  in¬ 
crease  the  angle  between  the  terrestrial  field  and  the  axis  of  the  magnet. 
For  angles  near  6  8  it',  the  solution  to  the  differential  equation  is 
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7 .  Procedure  and  Formulae  for  Evaluating  Engine  Performance-  Results. 

The  procedure  and  formulae  for  evaluating  the'  engine  perfor¬ 
mance  capability  are  summarized  in  this,  section.  The.  IKS.-  system,  of 
units  is  used  in  all  calculations . 


7.1  Satellite  Drag  (Equation  (5.3a)) 

If  the  satellite  is  assumed  to  have  a.  radius  a  and  a  length 
2a,  then  the  drag  is 

<Vv  *  (<Vv if  vv  <7-d 

Equation  (5.4)  gives  the  drag  coefficient' 

(O  »•  2.48. 

D  v 

We  take  as  satellite  dimensions  a  =  0.5  meters  ,  and;  satellite^ velocity 
vv  =  8000  raeters/sec.  Using  values  of  j>  from  the  ARDC:  model  atmosphere 
(sec  Appendix.  1)  the,  drag  force  on  the  vehicle  can  be  evaluated,  as 
a  function  of :  altitude.  The  results  are  shown  in  Figure  7.1.. 


7 .2  Ratio  ^ ^  of  Anode  Area  to  Solar  Cell  Area-  (Equation  5'.  13); 
We  have 


where 

*  *•  fraction  of  ambient  atmosphere  that  is  ionized, 

*  0.5 

v  *  thermal  velocity  of  electrons'. 

®  6' 

*  10  meters /sec 

I/Ij.  *  ratio  of  total  current  to  ion  current. 

=  2.0 

dP 

(tt)  3  solai  constant 
dA  s 

=*  1300  watts /sq^.raeter1 
•n  ■  solar  cell  efficiency 

CSC 

=  o.io 
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v  *  Alfven  speed  of  ions 
cr  4 

(v  >*1.6  x  10'  meters/ sec  (singly  ionized  Nitrogen) 

Cl»  ^ 

(v  )„*  2.4  x  10  meters/sec  (doubly  ionized  Nitrogen) 

-  cr  a. 

Using  the  values  of  y  from  the  ARDC  tables,  the  values  of  ^  can 
be  determined  as  function  of  altitude  for  both  values  of  the  Alfven 
speed.  Results  are  shown  in  Figure  7.2. 


7.3  Solar  Cell  Area  required  when  Cell  Array  Parallel  to  Line  of  Flight 
(Equation  5.8a) 


(FD>v 


3C  .  dP  .  T\ac 
1  dA;s  2 v 

cr 


1 


<v 


sc<  1  +*A> 


(7.3) 


where  from  Equation  (5.2a"> 


// 


<CD>sc 


2  u' 

T^S 

2  x  0.885 
1.772  x  7.9 


*  0.127 

The  other  parameters  are  given  above.  The  values  of  the  solar  cell  areas 
needed  for  the  two  Alfven  values  cf  the  exhaust  velocity  are  shown  as  a 
function  of  altitude  in  Figure  7.3.  This  calculation  indicates,  that  the 
system  will  not  operate  below  some  critical  altitude,  shown  here  to  be 
about  110  miles . 


7 .4  Anode  Area  Required  (Equation (5 . 14)) 


^f.A 

A  sc 


(7.4) 


Since  ^f.and  A  are  computed  above,  this  equation  can  be  evaluated 

A  SC 

immediately.  Values  of  the  anode  area  A^  are  shown  as  a  function  of 
altitude  in  Figure  7.4. 


7 • 5  Total  Power  Required  by  the  Propulsion  System  (Equation  5.5) 


P 


(  «  j 
'  HA -7 


dA’^s  IscAgc 


(7.5) 
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All  of  the  necessary  values  for  computing  P  are* above.  The  power 
requirements  are  shown,  as:  function  of  alitude  in  Figure  7;5. 

7.6  Total  Drag  on  the  Vehicle  (satellite  +  solar  cells  +  anode) 
(Equation  5.8) 


^FD*T  =  2v 


(7.6) 


Once  again  this  equation  can  be  evaluated  using  values  computed  above. 
The  results  are  shown  in  Figure  7i6. 

7.7  Total  Current  in  the  Electric  Discharge  (Equations  (5 .9), (2 . 2. 2),) 


where 


1*21, 

I  m.  m  v 

a  a,  cr 


Z  *  number  of  positive'  charges;  on  the  ion 

_ _ .  i  j _ »  _ _ _ '■*.  r  .  -,*“19  _ . _ 


(7.7) 


e  *  electronic  charge  *1.6-  x.  10*  coulombs 

•27 

nn  *  maos  of  the1  atom-  *14  x  A.  67  x  i0r  kg  for  Nitrogen, 

Using  the  values  for  the  total  drag  computed'  above  and  remembering  that 

Z  *  1  for  (v  ),,  and  Z  «  2  for  (v  )<»»  the  discharge  current  can  be. 
cr  i  cr  z  , 

computed  as  a  function  of  altitude.  The  results  are  shown  in  Figure  7.7. 

2‘ 

7.8  Minimum  Cross-Sectional  Area  for  the  Discharge  V R  in»  order  to  have 
adequate  mass  for  the  exhaust  beam.  (Equation  (5.7)-  ) 

(m)  -  e*p  v  tIR2  - 

cr  J  v  vcr 


so; that 


1  **Vt 


epv  v 
J  v  cr 


(7.8) 


or,  using  Equation  (7.7) 


•*  v 

a  .  I 


2h  lei 


(7.8a) 
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This  equation  can  be  immediately  evaluated  and  the  results  are 
shown  in  Figure  7.8.  It  should' be  emphasized  that  the  values  com¬ 
puted  are  the  minimum  values  for  the  cross-sectional  area  of  the 
discharge..  If  the  ionization  rate  is  such  that  less  than  1/2  of  the 
atoms  passing  through  the  discharge  are  ionized,  then  e<0.5  and  the 
cross-sectional  area  must  be  larger. 

7.9  Number  of  Turns  on  Magnet  to  produce  the  required  thrust . 

(Equation  (4.15)) 

*•  I*  (7-9) 

The  results  of  this  evaluation  are  shown  in  Figure  7.9.  The  numbers 
of  ampere- turns,  NI,  for  sthe  magnet  are  shown  in  Figure  7.9a. 


7.10  Mass  of  the  Magnet  (Equation  (4.18a).) 

M (7.10) 

w 

Once  the  material  for  the  magnet  has  been  selected,  this  equation  can 
be  evaluated.  Examination  of  Table  4.2  indicates  that  a  magnet  made  of 
sodium  clad  in  stainless  steel  would  likely  result  in  a  minimum  weight 
magnet.  For  this  case 

$  -  0.417  x)0~4  k«-fm 

°5  “ 

The  mass  of  the  magnet  is  shown  in  Figure  7.10  as  a  function  of  altitude. 

7.11  Surface  Area  of  the  Magnet  that  is  effective  in  radiation  cooling 
the  coil.  (Equation  (4.21b)) 

—  -  0.257  M^.  for  Sodium  (7.11) 

JW  J 

This  area  can  be  evaluated  using  the  magnet  mass  M  from  Section  7.10. 

The  magnet  surface  area  is  shown  in  Figure  7.11  vs.  altitude. 


■  <As>Min  ■  8( 
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7 .12  Magnet  Equilibrium,  Temperature  us Ing  the-  power,  dissipation  and 

solar  energy  as  inputs  and,  radiation  as  energy,  output,..  Assume,  that 

half  of  the  magnet  surface  is  directly  radiated  by  ther  sun., (Equation  (4.19)): 


P  .  ,  dP  .  (AsWn  /A  .  .J,. 

2  +  dA  2  ^AsWn  CR  & 


(7.12) 


where 


=  emissivity  of  wire  surface 

»  1.6 


o'  *  Stefan-Boltzmann  constant 
“  5.69  x,  10  ®  watts/meter^  °K^. 


The  magnet  surface  temperature  is  shown  in  Figure  7.12  vs.  altitude.. 


7.13  Volume  of  the-  Discharge  Necessary  to  produce  ionsatthe.  rate  re¬ 
quired  through  inelastic  collisions  of  electrons  and  atoms.,  (Equation  (2.11)) 


Vol‘  "  Z  |e|-nevepa  (oc/p)t 


C7..13) 


where 


Ij.  *  ion  current  in  exhaust  beam 


ng  *  average  electron,  density  in  volume  of  discharge 


p^  ■  ambient  pressure  of  atoms 


(ot/p)«  experimentally  determined  ionization  parameter  for  air 


*  7.6  (ion  pairs/mj  /  (newton/nr) 

vg  *■  electron,  thermal  velocity 

*  10**  meters/sec. 


Therefore  ye  have 


Vol.  «  164.5 


,  3 

meters  ., 


13 .  Discussion  ,of  the  Calculations 

Several  important  inferences  can  be  drawn  from  the  calculations . 
dne  of  ihe  most  important  is  that  the' use  of  solar  energy  to  power  the 
Advanced  Electric  Thruster  will  result  insufficient  thrust  for  drag 
make  up  at  altitudes  only  above  about  100  miles.  It  should  be  em¬ 
phasized,.  however,  that  the  engine  will  work  well  and  probably  better 
at,  lower  altitudes,  but  that  some:  source  of  power  other  than  or  in 
addition  to  solar  power  would'*  be  needed.  The  discharge  will  operate 
at  a  potential  of  between  25  and  30  volts  under  the  conditions  speci¬ 
fied  In  the  calculations .  This  matches  well  with  voltages  available 
from  solar  cell  arrays.  Further,  if  the  magnet  and  discharge  ars 
operated' in  series,  then  the  discharge  will  have  a  strong  positive 
characteristic  and  can  be  connected  directly  to  the  power  source, 
without  power  conditioning. 

The  performance  has  been  calculated  oh  the  assumption  that 
no  shock  wave  occurs  due  to  the  very  large  value  of  the  mean  free 
path  (more  than  1  meter)  bf'  the  gas  particles  in  the  altitude  range 
under  consideration.  In  fact,  effects  due  to  shock  wave  formation 
would  improve  engine  performance:  for  example,  increased  density  would 
improve  ionization  efficiency . 

The  maximum  axial  magnetic  field  strength  heeded  on  the  axis 
of  the  coil  is  about  57  gauss  for  accelerating  singly  charged  ions, 
and  about  43  gauss  for  accelerating  doubly  charged  ions.  These  field 
strengths  are  independent  of  the  altitude,  as  indicated  by  Figure  7.9a. 
In  practice,  there  is  little  or  nothing  that  can  be  done  to  restrict 
operation  to  one  of  either  the  singly  ionized  of  the  doubly  ionized 
mode.  As  the  ambient  density  decreases,  the  electron  temperature  will 
tend  to  rise,  leading  to  the  onset  of  higher  multiple  ionizations. 


9.  Pulsed  Operation 

The  Advanced  Electric  Thruster  can  be  operated  in  a  pulsed  mode 
as- well  as  steady-state. ,  However,  several 'problems  are  introduced  by 
pulsed  operation: 

1.  The  thrust  level  must  be  higher  by  the  ratio  of  the.  or¬ 
bit  .period  to,  the  operating  time  d.uring.  the, orbit.  This 
requires,  a  larger  cross-sectional  area  of  interaction  in 
order  to  intersect  enough  propellant,;  This  can  be  accom¬ 
plished  by  using  a  magnet,  with  a^ larger  magnetic  moment. 

2.  The  magnet  current  must  be  operated  in,  a  pulsed  inode  as 

well  as,  the  discharge  in  order  to  use  the  power  efficiently. 
.Ibis  can  be  accomplished  by  using  a  condenser  to  store  the 
magnet  power  and  arranging  to  open  the  circuit  between  the> 
magnet  coil  and- the  condenser. when,  all  of  the  energy  is 
stored  back  in  the  condenser;.  ,  „ 

There  is  a, potential  advantage  possible,  from  pulsed  operation,  if 
the  satellite  is  operating ,, in  the  ionosphere,  it  may  be  possible  that  a 
.considerable  number;  of  charged  particles  will;  be  accumulated  by.  the  mag-! 
netic  field  and  carried  along.  The  particles  will  be  accelerated  in  the 
exhaust  beam  when  the  discharge  current  .is  turned'  on. 
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10.  Integration  of  the  Advanced  Electric  Thruster  into  the  Spacecraft 

10 . 1  Possible  Configurations 

In  order  to  ensure. Mnimum  interference, -between- sensors,  etc.  on  the 
Satellite  and  the  AET  components,  the  thruster  will  be  positioned  on 
the  rear  of  the  satellite  or  if  desired  mounted  tip  to  0.5  meters  from 
the  body,  of  the  satellite.  The  coil  outer  diameter  should  be,. approx¬ 
imately  equal  to  the  diameter  of  the  satellite,  an^t  the  anode  should 
have  a  similar  . diameter.  An  .insulator  plate  placed  concentrically 
-around  the  cathode  will  extend  radially  out  to  the  anode  radius  .to 
prevent  metal  components  in  the  satellite  fronr  shorting  out  the  radial 
electric  field  between  the  electrodes.  As  mentioned  above;  an  axial 
gap-  . between  the  insulator  and  the  anode  ring  will  be;  necessary  to 
prevent  .the  flow  of  surface  currents  which  .would'  quickly  erode  the 
insulator. 

10.2  Mechanical  and  Electromagnetic  Interference  Effects 

Placing;  the  Advanced1 Electric  Thruster  at  or  beyond  the  rear  of 
the  satellite  should  make  possible  designing  out  any  mechanical  interfer¬ 
ence  problems  between  AET  and  satellite  components. 

The  d.c.  magnetic  field  of  the  engine  might  offer  some  problems 
of  interaction  with  moving  ferromagnetic  components  of' .the  satellite, 
and.  possibly,  with  electromagnetic  devices,  especially  under  conditions 
of  the  field  reversals.  There  are  various  ways  of  shielding  the  sus¬ 
ceptible  devices ferromagnetic  shielding,  suitable  design  and  alignment 
of  components,  and  so  on,,  or;  even  alternating  engine  operation  with 
that  of  the  device  or  devices... 

Alteration,  to  some  degree  of  the  engine  magnetic  field  by  the, 
presence  of  ferromagnetic  etc.  materials  in  the  vehicle  should  not 
degrade  engine  performance  appreciably,  and  may  possibly  be  so  designed 
as  to  enhance  performance .  ' 
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li.  Requirements  for  Wind-Tunnel  and/or  Flight  Test  of  the  Advanced 
Electric  Thruster  (AETV  Concept 

11,.  1  Critical  Engine  Parameters  Requiring  Investigation 

11.1.1  Engine  Components 

The  engine  consists  of  three  major  components:  the  anode-, 
the  magnet  coil,  and  the  cathode.  The  design  parameters  for  the  first 
two  are  well -known  and  can  ‘be  fabricated  once  the  material  for  con¬ 
strue  tionis  chosen.  However,  some  preliminary  experimental  work  should 
be  undertaken  to  investigate  cathode  design  parameters.  These  should 
determine  the  following: 

1.  The  best,  method  of  ensuring  that  the 'cathode  coating  material 
is  injected:  into 'the  cathode  cavity  at  the  proper  rate. 

2..  Which  coating  material  should  :be  used. 

3.  Therclationbetween,  the.mass.lossratefrom  the  cathode, and 
■the  current. 

4.  The  minimum  -power,  input  necessary,  to  .ensure  proper  cathode 
performance. 

5 .  Optimum  radiation  shielding  configurations . 

11, 2--' Performance  Parameters 
il.lj2.-l  Ionizing  Efficiency 

The  most  critical  phenomena  relative  to  the  engine  operation  that 
requires  experimental  investigation  is  the  effectiveness  with  which  the 
ambient  material  is  ionized  by  the  electric  discharge.  To  assess  this,  it 
is  necessary  to  conduct  a  test  in  an  environment  with  the  ambient  density 
close  to  that  found  at  the  altitudes  of  interest,  i.e.  those  over  100  miles. 

-Il.j-j2.2-  Anode  Operation 

Anode  shape,  size  and  positioning  require  experimental  study.  The 
anode  must  be  placed  in  such  a  position  as  to  prevent  the  discharge  from 
enveloping  the  vehicle,  i.e.  the  discharge  current  must  be  confined  to 
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the -space  behind1  the  plane  of  the  magnet.  The  outer  surface  of  the  anode 
ring  will  undoubtedly  have  to  be  covered  with  insulating  material  to 
ensure  that  this  does  occur.  Anode  shape  and  size,  and  other  factors 
involving  positioning  would  require  study  for  performance  analysis 
and  optimization. 

11.2  Scaling. 


There  are  a  number  of  non-dimensional  parameters  associated  with 
•the  performance  of  a  device  of  the  AET  type.  Some  are: 

1.  The  ratio  of  the  electron  and  ion  cyclotron  frequencies 
to  the  collision  frequency: 


U)  T  — - — 

e  6  me  nl^e 


and 
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(11.1) 

(11.2) 


2.  The  ratio  of  the  electron  and= ion  cyclotron  radii  to 
the  vehicle  radius,  a:  . 
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(11.4) 


3.  The  ratio  of  the  magnetic  pressures  to  the  gas  pressure: 

»2 


B 

For  the  applied  magnetic  field:  —z - 

2/oP 


(11.5) 


I.  -fie1* 


For  the  induced  magnetic  field:  p 

cr 

4.  The  ratio  of  the  ionizing  mean  free  path  to  the  vehicle 
radius : 


(11.6) 


n  q‘v  a 
a  ea  e 


(11.7) 
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5.  The  ratio  of  the  atom-atom  mean  free  path  to  the 
vehicle  radius: 


6.  The  Kach-  number: 


(11.8) 


(H.9) 


.7.  The  Reynold's  number: 

pva 

R  (11.10) 

e  ji  - 

Because  there  is  such  a  large  number  of  these  parameters^  design¬ 
ing  a  reduced  scale  experiment  would  be  very  difficult  ;  hence  it,  appears 
that  a  full  scale  experiment  would  be  the  only  feasible  method  of  testing 
the  concept. 


11.3  Size  of  the  Testing  Vacuum  Tank 

The  tank  should  be  sufficiently  large  that  the  magnetic  field 
will  drop  to  values  at  least  as  low  as  that,  of  the  earth's  magnetic 
field, at  the  wails.,  if  the  experiment  is  placed  at  the  center  of  the 
tank,  minimum  values  of  tank  length  and  diameter  can  be  termined  by' 
specifying  the  minimum  value  of  the  ampere  turns  of  the  magnet .  The 
calculations  in  Section  7  indicate  that  NI  should  be  about  400  ampere 
turns,  so  that  using  Equation (4.7)  ,  the  tank  length  2aZ  can  be  found: 


D  .=  ho'rtNI  .2 

e  ^IST  —  “  z2}3/2  (11.11) 

A 

where 

► 

Bg  **  earth's  magnetic  field 
rf  0.5  x  10"4:  Tesla 
NI  *  400  ampere  turns 
a  *  radius  of  test  equipment 
*  0.5  meter 
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Hence 


<1  +  3/2 


6j^^mjO0jlx400 
0.5  x  0.5  x  10-4 


126.2 


and 


V  4*92 

so  that,  since  2a=l,  the  tank  length  should  be  no  less  than  4.92  meters 
(16ft) . 

For  a  dipole  magnet,  the.  far  field  radially  out,  in  the  plane  of 
the  magnet  is  minus  oner  half  its  value  on  the  axis  for  the  same  distance 
from  the  center  of  the  coil.  .  The  tank  diameter  can  hence  be  estimated 
by  the  equation 


_  Uo'HNI  * 

Be  2  a  (1  +',RJi)3/k 


(11.12) 


or 


(1  +  R2)  3/2  =63.10 


R  »  3.85. 

Thus  the  tank  diameter  should  be  more  than  3.85  meters  (12.5ft)  (for 
a  vehicle  of  1  meter  diameter). 


11.4  The  Test  Environment 

The  test  engine  must  be  placed  in  an  air  flow  where  the  mass, 
momentum,  and  energy  fluxes  are  as  close  as  possible  to  that  encountered 
by  satellites  at  altitudes  between  iOO  and  300  miles .  An  electron  density 
of  between  10  and  10  electrons/  meter  needs  to  be  present  in  the 
flow. 

In  the  altitude  range  mentioned,  the  mean  free  path  varies  from 
5  meters  to  23,000  meters.  Because  of  this  and  since  the  size  of  vehicle 
under  consideration  is  1  meter,  little  or  no  significant  shock  wave  effects 
will  be  present,  .such  as  increases  in  density,  pressure  and  temperature 
over  the  ambient  values,  even  though  the  vehicle  would  be  travelling  at 
Mach  numbers  between  8  and  10. 


49 


Accurate  simulation  of  this  high-altitude  flow  field  is  not 
feasible  in  existing  facilities.  Figure  11 .1  illustrates  an  experimental 
configuration  which  should  offer  a  satisfactory  approximation  to  the 
conditions'  required.  In  it,  an  arc  heater  with  a  conical  nozzle  expands 
the  flow  to  a  Mach  number  of  between  2  to  4,  and  opens  to  a  large 
vacuum  tank  with  the  best  available  pumping  capacity,  so  that  the  flow 
will  con£S^ue  to  expand ,  decreasing  the  density  in  the  flow  to  as  low 

i 

a  value  as  [possible.  If  it  is  assumed  that  the  flow  crossr sectional 

j  "  . .  ’’  2 . 

area  at  the  test  position  is  about  18  meters- ,  then  the  mass  floy  rate 

of  air  through  the  arc  heater  should  yary  from  0.16  gm/sec  to  3.18  x  10  ^ 
gm/sec.  The  power  required  by  the  arc  heater  yould  be  frpo  about  20  kw 
down  ito  less  than  100  watts . 

It  appears  at  present  possible  only  to  produce -wip.d?  funnel  flows 
simulating. the  satellite  environment  ft  altitudes  under  10Q1  miles  dpe  to 
the  limitations  of  the  available  vacuum  facilities  @nd  tbs  technology  of 
producing  extremely  low  density  flows, 


11.5  .Test  Facilities  Available 

Because  existing  facilities  in  the  JJSA  do  npt  appear  capable  of  the 
accurate  highraltitude  flay  simulation  desired  for  the  full  scale  test, 
certain  compromises  and/or  auxifliary  devices  would  be  necessary,  The 
following  facilities  appear  to  be  the  best  available  for  consideration 
for  testing  the  Advanced  Electric  Thr^stef  . 

11,5.;1  A.E.D.C.  Low  Density  Tunnel  M 

This  tunnel  has  a  nitrogen  flay  from  gp  arc  heatpr  expanding  through 

I 

a  conical  nozzle  into  §  text  section  tank  3  ft  in  diamgtfr,  A  balance  is 
available  for  measuring  dffg  or  thrust  fppm  Q  to  100  milli-ppund§,  The 
tunnel  hp§  the  disadvantage  that  the  mass  flux  rate,  op,  is  more  than 
10^  too  high?  Furfhep,  the  tfnk  is  jtpp  small,  and  f ipcp  pumping  is 
accpmplish§4  by  ait  injegtprs  it  is  uplikely  that  f  loy  enough  density 
could  be  'achieved  fpr  meaningful  tests. 
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11.5.2  TA.E.D.C.  Aerospace  Chamber  C10V) 


This  chamber  is  designed  for — among  other  things — space  propulsion 
testing.  The  chamber  pressure  characteristics  (shown  in  Figure  11.1, 
as  taken  from  Reference  )  indicate  that  the  density  can  be  maintained 
at  a  level  of  10  torr  at  the  flow  rate  anticipated  from  the  arc  heater. 
This  is  about  1  order  of  magnitude  too  high  for  simulation  of  100  miles 
altitude.  The  tank  is  10  feet  in  diameter,  compared  with  the  minimum 
of  12.5  feet  desired.  If  under  certain  compromises  this  tank  were  used, 
a  small  nitrogen  arc  heater  would  have  to  be  designed  and  built  as  part 
of  the  experiment. 

11.5.3  A.E.D.C.  Aerospace  Environmental  Chamber  (Mark  I) 

This  tank  is  approximately  34  feet  in  diameter  and  65.5  feet  long, 
and  so  is  clearly  large  enough.  It  also  has  a  pumping  capacity  adequate  to 

-7 

conduct  a  test  at  a  pressure  of  10  torr,  simulating  an  altitude  of 
100  miles.  The  maximum  throughput  of  nitrogen  at  this  pressure  would  be 
only  a  few  milligrams/sec.  This  means  that  the  arc  heater  must  operate 
at  as  low  a  power  as  feasible,  probably  about  1  kw.  Once  again,  if  this 
tank  were  used,  the  arc  heater  and  balance  would  have  to  be  designed  and 
built  as  part  of  the  experiment. 

11.5.4  N.A.S.A.  Lewis  Electric  Propulsion  Facility 

Meaningful  tests  could  be  conducted  in  this  facility  if  it  could 
be  made  available.  It  is  possible  that  one  of  the. NASA  Lewis  MPD  arcs 
could  be  used  as  the  arc  heater.  Balances  are  available  to  measure  the 
drag  and  thrust. 

11.5.5  Other  Facilities 

There  are  a  number  of  facilities  in  industry  and  at  other  govern¬ 
ment  installations  that  might  be  considered,  such  as  the  E0S-USAF 
electric  propulsion  facility,  and  the  WPAFB  electric  propulsion  facility. 
However,  they  appear  well  below  the  previously  mentioned  facilities 
in  meeting  standards  of  suitability  for  the  experiment  as  described. 
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12.  Proposed  Experimental  Program 

The  following  program  is  reconmended  as  the  next  phase  of 
this  study.  It  is  so  constructed  that  each  section  could  be  done 
as  an  independent  study. 

12.1  Cathode  Studies 

A  cathode  structure,  similar  to  that  shown  in  Figure  3.3,  should 
be  designed  and  built,  and  tests  conducted  in  a  small  vacuum  chamber  to 
determine  relations  among  the  following  variables: 

1.  temperature  of  the  cathode  cavity 

2.  material  used  to  coat  the  cathode  . 

3.  mass  loss  rate  through  the  orifice 

4.  strength  of  the  applied  solenoidal  magnetic  field 
all  as  functions  of  the  total  current  carried  between  the  cathode  cavity 
and  an  anode  placed  at  least  10  cm  from  the  cathode. 

12.2  Investigation  of  Engine  Performance  in  a  Flow  Field  with  a  Low  Level 
of  Ionization 


This  experiment  could  be  conducted  on  a  transient  pressure  basis, 
in  that  the  tests  would  be  initiated  in  the  tank  at  its  highest  vacuum  and 
data  accumulated  before  jet  influxes  increased  pressure  excessively.  The 
objective  would  be  to  determine  the  extent  to  which  the  accelerator  can 
utilize  the  ionized  fraction  of  gas  at  the  atmospheric  E- layer  as  propellant 
and  to  provide  discharge  breakdown  nucleii.  The  tests  would  also  determine 
variation  of  the  discharge  as  ambient  pressure  is  reduced. 

The  facilities  and  components  needed  to  conduct  this  series  of 
tests  would  be 

1.  a  vacuum  chamber  at  least  10  ft  in  diameter  that  has 
a  pumping  system-capable  of  reducing  the  tank  pressure 
to  under  10  torr. 

2.  a  small  (l-10kw)  axfc  jet  that  can  operate  on  a  pulsed 
basis  (O-lOsec).  This  would  be  used  as  a  pulsed,  point 
source  of  gas  and  plasma.  The  conical  expansion  of  this 


flow  will  reduce  electron  density  and  mass  flux  at  the 
test  engine  to  the  desired  values. 

3.  a  test  propulsion  engine  consisting  of  a  magnet  coil) 
an  anode  assembly,  a  cathode  assembly  and  a  battery 
power  supply  of  approximately  1  kw  capability. 

These  tests  would  also  determine  the  feasibility  of  operating  the  engine 
by  acceleration  only  of  the  existing  ions  in  the  flow,  and  determine  the 
current  and  voltage  levels  at  which  this  can  be  done,  as  a  function  of  the- 
magnetic  field  strength  and  the  flow  parameters. 

12.3  Investigation  of  Engine  Performance  in  Simulated  Flow  Field 

This  series  of  tests  would  place  the  thruster  on  a  balance  to 
measure  drag,  thrust,  and  torque.  The  tests  would  be  conducted  in  a 
facility  that  could  simulate  the  required  mass  flux  rate  and  electron  on 
a  steady  basis  (i.e.  for  1  minute  or  more). 
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Figure  2.1  Schematic  of  Advanced  Electric  Thruster  (SERJ) 
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Fifl.12.4  Aerospace  Chamber  (10V)  Operating  as  a  Space  Propulsion  Test  Facility: 
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TABLE  4.2 


Material 

t,  , 

kgm/mJ 

1/cr 

ohm-m 

p/o' 

kgm-ohm/m 

Silver 

10.5  xlO3 

1.63xl0*8 

1,7  : 

«10*4 

Copper 

8.89  " 

1.72 

fl 

1.53 

II 

Aluminum 

2.70  ” 

2.69 

It 

0.725 

II 

Magnesium 

1.74  " 

4.46 

ft 

6.778 

II 

Lithium 

0.534  " 

8.55 

ft 

0.457 

II 

Sodium 

0.971  " 

4.30 

tl 

0.417 

II 

Potassium 

0.870  " 

6.10 

II 

0.531 

II 
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Relevant  parameters  of  the  ARDC  STANDARD  A1H0SPHERE 
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Table  1-4  (coat.) 
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SCALE 
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DENSITY  j 

FREQ. 

FREE  PATH 

Z;m 

H,m' 

km 

HBSB 

EEC9 

BE3H 

L,m 

B.v  i  A'  '"9HHI 

63362 

7.0484 

4.5546+  21 

415.36 

1.1196+  6 

3.7093  -  4 

60000 

65322 

6.7893 

3.5435 

407.53 

8.5476  +  5 

4.7677 

68000 

67280 

2.7301 

■  399.55 

6.4564 

6.1883 

70000 

69238 

6.2706 

2.0813 

391.41 

4.8217 

8.1175 

72000 

71194 

1.5685 

383.10 

3.5567 

1.0771-3 

74000 

73148 

5.7512 

:  1.1674 

374.61 

2.5866 

1.4472 

76000 

75102 

5.491 

8.571  +.20 

365.9 

1.857 

1.971 

78000 

77054= 

5.231 

6.199 

357.1 

1.310 

2.725 

80000 

79006 

4.972 

4.410 

348.0 

9.062  +  4 

3.831 

82000 

80956 

4.975 

2.849 

348.0 

6.075 

5.728 

84000 

82904 

4.978 

1.973  +  20 

348.0 

4-064  +  4 

8.562  -3 

84852 

4.981 

1.321 

348.0 

2.730 

1.279  -2 

86798 

4.984 

8.840  +  19 

348.0 

1.821 

1.911 

88743 

4.987 

5.918 

348.0 

1.219 

2.855 

32000 

90687 

5.073 

3.904 

350.9 

8.108  +  3 

4.327 

94000 

92630 

5.310 

2.540 

<  358.9 

5.395 

6.652 

96000 

94572 

5.548 

1.684 

366.7 

3.655 

fE* 

98000 

96512 

5.786 

1.136 

374.3 

2.516 

1.488 

98451 

6.023 

7.783  +  18 

381.8 

1.759 

2.171 

102000 

100389 

6.261 

.  5.413 

389.2 

1.347 

3.121 

104000 

102326 

6.500 

3.816  +18 

396.4 

8.953  +  2 

4.428  -1 

104261 

6.738 

2.724 

403.5 

6.506 

6.202 

108000 

106196 

7.555 

l;830 

427.1 

4.636 

9.233 

108129 

8.731 

1.240 

459.0  • 

3.368 

1.363  +  0 

110061 

9.908 

8.823  ♦  17 

488.8 

2.553 

1.915 

114000 

111991 

11.09 

6;  525 

516.9 

1.996 

2.589 

116000 

113921 

12.26 

4.975 

543.5 

1.600 

31396 

118000 

115849 

13.44 

3.890 

568.8 

1.310 

4.343 

117777 

14.62 

,  3.105 

593.1 

1.090 

5.441. 

125000 

122589 

17.57 

1.899 

649.7  ; 

7.301  ♦  1 

8.898 

130000 

127395 

20.53 

1.254  +  17 

701.7 

5.306  +  1  ■ 

1.347  ♦  1 

135000 

132192 

23.49 

8.764  +  16 

750.6 

3.891 

1.928 

140000 

136983 

26.46 

,  6.395 

795.3 

3.010 

2.642 

145000 

141786 

29.43 

4.829 

838.1 

2.395 

3.499 

146542 

32.40 

3.748 

878.8 

1.950 

155000 

151310 

35.38 

2.977 

917.6 

1.617 

5.675 

160000 

156071 

38.36 

2.411 

954.7 

1.363 

7.007 

165000 

160825 

1.998 

987.4 

1.168 

8.456 

170000 

165572 

42.62 

1.722 

UA5. 

1.034 

9.813 

175000 

170311 

44.11 

1.495 

1021. 

1.130  +  2 

175043 

44.91 

1.324  +  16 

1030, 

<  8.070  ♦  0 

1.276  +  2 

135000 

179768 

45.71 

1.177 

w# 

7.234 

1.435 

190000 

184485 

46.51 

1.050 

1046. 

6.503 

ft 

195000 

189196 

47.31 

9.385  +  13 

1055. 

5.858 

200000 

193899 

48.12 

8.406 

1063. 

5.288 

2.010 

198595 

48.92 

7.545 

1071. 

4.762 

2.239 

210000 

49.58 

6.805 

1077, 

4.339 

2.483 

.  50. 16 

6.154 

1083. 

*  3.944 

2.745 

212640 

50.75 

5.573 

1068. 

3.589 

3.032 

225000 

217307 

51.34 

5.052 

1094. 

3.270 

:  3.344 

¥ 


£ 

I' 


* 

s 


1 

i 


t 


i 


b 


¥ 


h 

> 

l 


j 

1 


\ 

i 

> 


i 

) 

1 

I 


MODEL  ATMOSPHERES 

Table  1-2  (cont.) 


ALTITUDE 

TEMPERATURE 

PRESSURE 

- r  "T 

Z,m 

H,m‘ 

T,°K 

■SSI 

P,mb 

HH 

P,mm  Hg 

64000 

66000 

68000 

70000 

72000 

74000 

76000 

78000 

80000 

82000 

63362 

65322 

67280 

69238 

71194 

73148 

75102 

77054 

79006 

80956 

236.03 

227  .  21 

218. 40 

209. 59 

200. 79 
191.99 

183.2 

174.4 

165.7 

165.7 

236.03 

227.21 

218.40 
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200.79 
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174.4 

165.7 
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1.4838  -  1 
1.1113 

8. 2298  -  2 
6.0209 
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3.0937 
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1.492 
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1.5131  +  0 
1.1332 
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6.1396 

4.4327 

3.1547 
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1.522 
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1.1130  -  1 
8.3354  -  2 
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4. 5160 

3.2605 
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1.626 

1.119 

7.563  -  3 

5.058 

84000 

86000 
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90000 
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96000 
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102000 

82904 

84852 

86798 

88743 

90687 

92630 

94572 

96512 
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100389 
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168.4 
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6.806  -4 
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1.158 
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108129 
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111991 

113921 

115849 
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325.0 

363.1 
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439.1 
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249.6 
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326.9 

365.5 

404.1 

442.6 

481.2 
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8.341-5 

6.274 
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2.358 
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6.395 
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3.614 

3.207 

2.858 

8.632i-  6 

6.876 

5.627 

4.704 

4.001 

3.452 

3.014 

2. 658 
2.359 

2. 102 

180000 

185000 

190000 

195000 

200000 

205000 

210000 

215000 

220000 

225000 

175043 

179768 

184485 

189196 

193899 

198595 

203284 

207965 

212640 

217307 

1371. 

1381. 

1389. 

1397. 

1404. 

1411. 

1414. 

1414. 

1414. 

1414. 

1451. 

1474. 

1498. 

1522. 

1545. 

1569. 

1587. 

1604. 

1620. 

1636. 

2.505  -8 
2.243 

2.013 

1.809 

1.622 

1.470 

1.328 

1.201 

1.088 

9.863  -  7 

2.554  -5 

2. 287 

2.052 

1.845 

1.661 

1.499 

1.354 

1.225 

1.109 

1.006 

1.879  -  6 
1.682 
1.510 
1.357 
1.222 
1.102 
9.959  -  7 
0.009 
8.159 
7.398 
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Table  1-2  (cont.) 


ALTITUDE 

TEMPERATURE 

PRESSURE 

Z,m 

T,°K 

tm-  °k 

P,mb 

P,kgf  m~2 

P.mm  H 

230000 

221968 

1415. 

1653. 

8. 953  -  7 

9. 129  -  6 

6.715-7 

240000 

231267 

1415, 

1685. 

7.401 

7.547 

5.551 

250000 

240539 

1415. 

1718. 

6.143 

6.265 

4.608 

260000 

249782 

1416. 

1750. 

5.120 

5.221 

3.841 

270000 

258998 

1417. 

1782. 

4.  284 

4.369 

3.213 

280000 

268185 

1418. 

1814. 

3.598 

3.669  i 

2. 699 

290000 

277345 

1420. 

1846, 

3.033 

3.092 

2,  275 

300000 

286478 

1423. 

1878, 

2.  565 

2.615 

1.924 

310000 

295583 

1426. 

1910. 

2.176 

2.219 

1.632 

320000 

304661 

1430. 

1942, 

1.853 

1.889 

1.390 

330000 

313711 

1435. 

1974. 

1.582  -  7 

1.613  -  6 

1.187-7 

340000 

322735 

1440. 

2005. 

1.355 

1.382 

1.016 

350000 

331731 

1445. 

20.J7. 

1.164 

1.187 

8.729  -8 

360000 

340701 

1451. 

2068. 

1,002 

1.022 

7.518 

370000 

349644 

1458. 

2099. 

8.656  -  8 

8.827  -  7 

6.492 

380000 

358561 

1465. 

2131. 

7.495, 

7.643 

5.622 

390000 

367451 

1473. 

2162. 

6.506 

6.634 

4. 880 

400000 

376315 

1480. 

2193. 

5.661 

5.773 

4;  246 

410000 

385152 

1489. 

2224. 

4.938 

5.035 

3.703 

420000 

393964 

1497. 

2255. 

4.316 

4.401 

3.238 

430000 

402749 

1506. 

2285. 

3.782  -  8 

3.856  -  7 

2.837  -6 

440000 

411509 

1516. 

2316. 

3.320 

3.386 

2.491 

450000 

420243 

1525. 

2347. 

2.921 

2.979 

2.191 

460000 

428951 

1535. 

2377. 

2.576 

2.626 

1.932 

470000 

437634 

1545. 

2407. 

2.  275 

2.320 

1.707 

480000 

446291 

1555. 

2438. 

2.014 

2.053 

1.510 

490000 

454923 

1566. 

2468. 

1.786 

1.821 

1.339 

500000 

463530 

1576. 

2498. 

1.586  . 

1.617 

1.190 

510000 

472111 

1587. 

2528. 

1.412 

1.439 

1.059 

520000 

480668 

1598. 

2558. 

1.258 

1.283 

9. 438  - S 

530000 

489200 

1609. 

2588. 

1.123  -  8 

1.146  -  7 

8.427  -4 

540000 

497707 

1621. 

2618. 

1.005 

1.025 

7.536 

550000 

506189 

1632. 

2647. 

9. 000  -  9 

9.178  -  8 

6.751 

560000 

514647 

1644. 

2677. 

8.075 

8.234 

6.056 

570000 

1655. 

2706. 

7.255 

7.398 

5.442 

580Qfm 

51*489 

1667. 

2736. 

6.528 

6.657 

4.896 

5900C. 

539874 

1679. 

2765. 

5.882 

5.998 

4.412 

600000 

548235 

1691. 

2794. 

5.308 

5.413 

3.981 

610000 

556571 

1703. 

2824. 

4.796 

4.891 

3.597 

620000 

564884 

1715. 

2853. 

4.339 

4.425 

3.255 

630000 

573173 

1727. 

2882. 

3.931  -  9 

4.009  -  8 

2.949  -1 

640000 

581438 

1739. 

2911. 

3.566 

3.636 

2.675 

650000 

589680 

1751. 

2940. 

3.239 

3.303 

429 

660000 

597898 

1763. 

2968. 

2.945 

3.003 

2.209 

670000 

606092 

1775. 

2997. 

2.681 

2.734 

2.011 

680000 

614263 

1788. 

3026. 

2.444 

2.492 

1.833 

690000 

622411 

1800. 

3054. 

2.  230 

2.274 

1.673 

700000 

630536 

1812. 

3083. 

2.037 

2.077 

1.528 

V 
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Table  1-3  (Cont) 


ACCELERATION 

SPECIFIC 

MOLECUL 

ALTITUDE 

OF  GRAVITY 

WEIGHT 

DENSITY 

WEIGHT 

Z;  m 

H,  m 

g,  m  sec*^ 

u,  kgf  m"^ 

P,  kg  m'3 

M 

230000 

221968 

9.134 

1.758  -  10 

1.887  -  10 

24;.?9. 

240000 

231267 

9. 100 

1.421 

1;  530 

24.32 

250000 

240539 

'9.078 

1.154 

1.246 

23.87 

260000 

249782 

9.051 

9.408  -  11 

1.019 

23.44 

270000 

258998 

9.024 

7.706 

8.375  -  11 

23.03 

280000 

268185 

8.996 

6.338 

6,909 

22.65 

290000 

277345 

8.969 

5.233 

5.722 

22.28 

300000 

286478 

8.942 

4.338 

4.757 

21.95 

310000 

295583 

8.916 

3.609 

3.969 

21.63 

320000 

304661 

8.889 

3.013 

3.324 

21.33 

330000 

313711 

8.862 

2.524  •  11 

2.792  -  11 

21.06 

340000 

322735 

8836 

2.121 

2.354 

20.80 

350000 

331731 

8,810 

1.788 

1.991 

20.55 

360000 

340701 

8.783 

1.512 

T.  688 

20.33 

370000 

349644 

8.757 

1.283 

1.436 

20.12 

380000 

358561 

8.731 

1.091 

1.225 

19.92 

390000 

367451 

8.705 

9. 307  -  12 

1.048 

19.73 

400000 

376315 

8.680 

7.960 

8.994  -  12 

19.56 

410000 

385152 

8.854 

6.826 

7.736 

19.39 

420000 

393964 

8.628 

5.868 

6.670 

19.24 

430000 

402749 

8.603 

5.057  •  12 

5.765  -  12 

19.09 

440000 

411509 

8.578 

4.369 

4.995 

18.96 

450000 

420243 

8.552 

3.783 

4.338 

18.83 

460000 

428951 

8.527 

3.283 

3.775 

18.71 

470000 

437634 

8.502 

2.855 

3.293 

18.59 

480000 

446291 

8. 477 

2.488 

2.878 

18.48 

490000 

454923 

8.453 

2.173 

2.521 

18.38 

500000 

463530 

8.428 

1.901 

2. 212 

18.28 

510000 

472111 

6,403 

1.667 

1.945 

18.19 

520000 

480668 

8,379 

1.464 

1.714 

18.10 

530000 

489200 

8.355 

1.289  -  12 

1.512  -  12 

18.01 

540000 

497707 

8,330 

1.136 

1.337 

17.93 

550000 

506189 

8.306 

1.003 

1.184 

17.86 

560000 

514647 

8.282 

8.875  -  13 

1.051 

17.78 

570000 

523080 

8.258 

7.864 

9.339  -  13 

17.72 

580000 

531489 

8.235 

6.980 

8.313 

17.65 

590000 

539874 

8.  211 

6.205 

7.411 

17.58 

600000 

548235 

8.187 

5.525 

6.617 

17.52 

610000 

556571 

8.164 

4.926 

5.917 

17.47 

620000 

564884 

8.140 

4. 399 

5.299 

17.41 

630000 

573173 

8.117 

3.934  -  13 

4.753  -  13 

17.36 

640000 

581438 

8.094 

3.523 

4.268 

17,30 

650000 

589680 

8.071 

3. 159 

3.839 

17.25 

660000 

597898 

8.048 

2.837 

3.457 

17.21 

670000 

606092 

8.025 

2;  551 

3.117 

17.16 

680000 

614263 

8.002 

2.296 

2.814 

17.12 

690000 

622411 

7.979 

2.070 

2.544 

17.07 

700000 

630536 

7.957 

1.868' 

2.302 

17,03 

A  1-7 
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2.  PHYSICAL  PROCESSES 

Ail  electric  propulsion  engines  rely  upon  two  basic  processes: 

1.  Ion  production 

2.  Ion  acceleration 

In  accomplishing  these  ends,  certain  penalties  must  be  paid,  which  are 
usually  evaluated  in  terms  of  propellant  mass  utilization  and  electric 
power  utilization.  The  axisymmetric  Hall  current  accelerator,  or  MPD 
arc, accomplishes  all  of  the  three  processes  mentioned  above  by  means 
of  a  dc  arc  struck  between  electrodes  placed  symmetrically  in  a  sole- 
noidal  magnetic  field.  However,  the  measurements  that  can  be  made  are 
such  that  it  is  virtually  impossible  to  evaluate  propellant  mass  uti¬ 
lization  and  power  utilization  independently.  Thus,  great  simplicity 
of  design  is  gained  at  the  expense  of  great  complexity  in  the  inter¬ 
action  of  the  engine  operating  mechanisms. 

One  of  the  features  of  the  MPD  arc  that  appears  most  baffling  to 
physicists  is  the  apparently  uncorrelated  behavior  of  the  arc  current 
and  the  ion  flux  rate.  Since  the  electrons  can  carry  an  appreciable 
fraction  of  the  arc  current,  considerable  energy  from  the  electric 
field  is  transferred  initially  to  the  electron  internal  and  kinetic 
energy,  from  where  it  can  later  be  transferred  to  the  ions.  This 
allows  the  ions,  under  some  circumstances,  to  leave  the  engine  with 
energies  greater  than  they  could  achieve  by  falling  through  the  poten- 
fial  drop  of  the  discharge.  The  detailed  process  by  which  this  can  be 
accomplished  will  be  discussed  in  the  following  paragraphs. 
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2 . 1  Ion  Production 

The  momentum  conservation  equations  indicate  that  some  torque 
and  force  reaction  must  occur  on  the  engine  (see  Section  3).  For  this 
to  happen,  ions  must  be  produced  in  and  expelled  from  the  motor  with 
axial  and  angular  velocity.  The  fewer  the  ions,  the  higher  the  result¬ 
ant  velocities  must  be  and,  consequently,  the  beam  power  must  be  higher. 
On  the  other  hand,  as  the  ion  production  and  expulsion  rate  increases, 
the  power  in  ion  production  becomes  very  high,  hence  some  ion  flow  rate 
must  exist  at  which  the  electrical  power  into  the  discharge  is  a  mini¬ 
mum.  If  the  arc  current  is  Held  fixed,  this  implies  that  the  discharge 
seeks  a  minimum  voltage  mode  in  which  to  operate.  The  above  argument 
then  indicates  that  the  arc  accomplishes  this  end  by  ionizing  the  opti¬ 
mum  amount  of  propellant  to  expell  in  the  exhaust  beam.  The  key,  then, 
to  explaining  the  operating  characteristics  of  the  engine  lies  in 
determining  the  optimum  ion  exhaust  rate  and  in  understanding  the  pro¬ 
duction  process  for  these  ions  throughout  the  volume  of  the  discharge. 

Previously  (Ref.  1),  it  had  been  concluded  that  volume  ion¬ 
ization  away  from  the  electrode  surfaces  could  not  account  for  all  of 
the  ions  produced.  At  that  time,  it  was  assumed  that  the  electron 
temperature  was  about  10,000°K.  Since  then,  data  on  anode  heating 
have  indicated  that  the  electron  temperature  may  be  much  higher,  e.g., 
sa  30,000  -  50,000°K.  This  raises  the  probability  of  volume  ion  pro¬ 
duction,  through  electron  atom  collisions,  to  values  where  it  can 
account  for  most,  if  not  all,  of  the  ions  used  in  the  exhaust  beam. 

The  atoms  are  not  confined  to  any  great  extent  by  the  elec¬ 
tric  discharge.  Hence  the  flow  field  of  the  gas  will  be  substantially 
the  same  as  one  would  find  for  the  gas  issuing  from  the  orifice  with 
no  discharge  present.  The  discharge  must  now  encompass  this  gas  and 
adjust  the  electron  temperature  and  density  throughout  its  volume  so 
as  to  ionize  the  optimum  amount  of  propellant.  Clearly,  to  get  good 
propellant  utilization,  the  injected  flow  rate  should  be  close  to  the 
ion  flow  rate  in  the  beam.  However,  in  the  interest  of  obtaining  best 
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overall  efficiency,  it  may  be  necessary  to  inject  slightly  more  pro¬ 
pellant  than  is  used  by  the  beam.  This  would  help  to  restrict  the 
volume  of  the  discharge  and  perhaps  keep  the  electron  temperature  to 
lower  values,  thus  reducing  power  loss  by  electron  energy  convection 
to  the  anode. 

Since  a  vast  majority  of  atom-electron  collisions  are  elastic 
collisions,  the  energy  used  to  ionize  one  atom  must  be  considerably 
greater  than  the  ionization  potential  of  the  atom.  The  energy  differ¬ 
ence  i-a  transferred  into  the  internal  energy  of  the  atoms  and  ions. 

If  most  of  the  injected  mass  is  eventually  ionized,  this  energy  is  not 
lost  but  is  available  to  be  transferred  into  beam  kinetic  energy  by 
eventual  expansion  througu  the  magnetic  nozzle.  It  is  obvious,  of 
course,  that  this  internal  energy  of  the  heavy  particles  can  never  be 
higher  than  that  of  the  electrons.  This  argument  indicates  that  the 
ionization  process  need  not  be  efficient.  However,  the  number  of 
inelastic  collisions  tnat  excite  the  electrons  to  states  that  can 
radiate  should  be  minimized.  This  is  basically  a  problem  of  propellant 
selection. 

The  electron  internal  energy  results  from  the  electron  current 
passing  through  the  potential  drop  and  being  randomized  by  collisions 
with  heavy  particles.  For  this  reason  it  would  be  expected  that  the 
highest  electron  energy  would  be  found  in  the  anode  sheath,  after  the  4 
electrons  have  fallen  through  most  of  the  potential  drop.  This  is 
fortuitous,  since  it. is  precisely  in  this  region  that  the  highest  pro¬ 
duction  rate  is  wanted,  to  let  the  ions  gain  a  maximum  of  kinetic 
energy  by  falling  through  the  potential  back  toward  the  cathode  jet 

2.2  Ion  Acceleration 

Momentum  can  be  transferred  to  the  ions  from  electric  fields 
or  from  collisions  with  other  particles.  For  convenience,  the  momentum 
exchange  processes  in  a  fully  ionized  gas  through  which  an  electric  dis¬ 
charge  is  passing  shall  be  discussed. 
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Locally,  the  force  ot;.  each  Ion  is  given  by  the  following 

J  X 


expressions: 

Axial: 

i  °  i  {e; 

Radial: 

Mfc 

Azimuthal: 

1  e  |  {w. 

1  e  I  -{e.  v_  B  —  w_  B.  - 

1  1  l  r  I  z  I  0 


9  c 


T  B  -  uT  B - 

I  r  I  z  a 


If  an  electric  discharge  is  established  in  a  uniform  axial 
magnetic  field  we  shall  call  the  region  where  the  current  flows  down¬ 
stream  the  anode  sheath  and  the  region  where  it  flows  upstream  the 

cathode  jet.  If  a  cathode  of  maximum  diameter  R  is  surrounded  by  an 

c  # 

anode  ring  of  diameter  R^  (R^  >  Rc)  we.  ask  the  questions: 

1.  Does  the  cathode  jet  expand  out  to  meet  the  anode  sheath? 

2.  Does  the  anode  sheath  contract  in  diameter  to  meet  the 
cathode  jet? 

3.  Do  both  (1)  and  (2)  occur  simultaneously? 

Consider  first  the  cathode  jet.  If  the  cathode  jet  is  to 
expand  outward,  conservation  of  momentum  states  that  the  axial  momentum 
of  the  jet  must  increase  and  the  rotational  momentum  of  the  jet  must 
increase.  However,  the  local  (E  +  v  x  B)  axial  and  tangential  electric 
fields  are  both  in  the  wrong  direction  to  accelerate  the  ions  and  the 
momentum  must  hence  be  transferred  to  them  by  electron  collisions. 

This  is  a  highly  dissipative  process,  resulting  in  strong  heating  of 
the  electrons.  This  increases  the  rates  of  entropy  production  over 
that  caused  by  ion-electron  drag  in  a  purely  dissipative  plasma  (no 
body  force) . 

In  the  anode  jet  the  situation  is  quite  different.  Both 
Ez  and  u^  Bg  are  in  the  positive  z  direction,  thus  helping  to  accel¬ 
erate  the  ions  axially.  The  only  dissipative  or  entropy  producing 
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i 
* 

!  1.  A  uniform  diameter  cathode  jet  is  established  which  carries 

only  a  small  fraction  of  the  injected  mass  flow  rate* 

2.  An  annulus  of  plasma  is  established  off  of  the  anode  face. 

The  injected  mass  is  accumulated  within  this  annulus.  'All 
of  the  discharge  current  passes  through  the  annulus. 

3.  The  average  radius  of  this  anode  sheath  decrease's  downstream. 

.This  causes  the  ions  in  the  sheath  to  be  accelerated 

slightly  in  the  axial  direction  and  to  be  spun'  up  to  high' 
azimuthal  velocities.  The  electrons  are  simultaneously 
heated,  mainly  by  the  ion-electron  drag  in  the’ azimuthal 
direction,  where  the  azimuthal  electron' motion -is  helping 
to  spin  up  the  ions.  " 

4i  The  anode  sheath  eventually  meets  ihe  cathode  jet  at  z  =*  L 
i  and  the  current  path  is  completed.  "No  discharge  current 

flowsat  axial  positions  beyond  this  point,- 
5.  At  positions  of  z  >  L,  the  magnetic  field  acts  Tike  a 

magnetic  nozzle.  As  the  field  diverges  tne  ions  are  accel¬ 
erated  axially  by  two  processes: 

a.  Conversion  of  angular  momentum  requires  that  as  the 
jet  radius  increases,  rotational  ion  energy  must  be 
1  transferred  to  axial  and-  radial  kinetic  energy. 

r 

;  b.  The  high  energy  electrons  tend  to  expand  out  of  the 

\ 

nozzle  ahead  of  the  ions,  thus  setting  up  a  positive 
axial  electric  field  that  accelerates  the  ions.  In 
this  manner,  all  of  the  energy  of  the  particles  in  the 
beam  can  be  converted  into  the  kinetic  energy  of  the 
ions.  Obviously,  some  of  tni.s  energy  will  reside  in 
the  radial  motion  with  the  result  that  the  expansion 
will  not  be  100  percent  efficient, 

i 

I  An  attempt  to  analyze  the  "heating"  (constant  magnetic  field) 

!  region  of  such  a  device  is  made  in  Section  4  of  this  report. 
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3.  INTEGRALS  OF  THE  MOMENTUM  EQUATIONS 

There  are  a  few  cases  where  the  net  electromagnetic  force  in  an 
axisymmetric  body  can  be  computed  without  detailed  knowledge  of  the 
distributions  of  current  and  magnetic  field.  In  the  general  case,  the 
momentum  equations  must  be  integrated  simultaneously  with  continuity, 
energy,  Ohm's  laws.  Maxwell's  equations  and  the  equations  of  state. 

The  cases  chosen  here  are  such  that  the  integrand  (force  per  unit 
volume)  can  be  put  into  the  form  of  a  divergence,  by  using  Maxwell's 
equations  and  simplified  momentum  equations.  These  forces  can  be 
integrated  in  terms  of  total  current,  radius,  and  applied  magnetic 
field. 

3.1  Pressure  Due  to  T  BQ  Pinch 

_ Z  o _ 

Here  the  average  pressure  on  the  cathode  is  computed  in 
*- 

terms  of  the  current  and  radius  of  attachment.  The  equations  used- 
are  a  momentum  equation 

5?  -  -  Ve-  Bfc'R)  -  p0  ,  U> 

an  induction  equation 

Vr-0>-°  •  <2> 

a  total  current  integral 

R 

I  -J  Jz(r)  2 nr  dr  ,  ■  (3) 

o 
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and  a  definition  of  average  pressure 


p  =  f  p(r)  2  nr  dr  (4) 

n  R  J 
o 

Combining  the  above  relations,  it  follows  that  independent 

of  the  distribution  of  J  ,  tht  average  cathode  pressure  is  given  by: 

2 


+  AJL 

Pav  “  ?0  8n2  R2 


(5) 


The  pressure  given  by  Eq.  5  will  act  on  the  cathode  to  give 

2 

a  thrust.  This  thrust  force  is  given  by  hq1  /8n  and  is  independent 
of  the  distribution  of  the  current  densif"*  at  the  cathode,  and  of  the 
size  of  the  cathode  attachment. 

3.2  Thrust  Due  to  Pumping 

The  amount  of  thrust  in  an  axially  symmetric  volume  due  to 
radial  currents  and  induced  azimuthal  magnetic  fieici  can  v  i  evaluated 
in  terms  of  the  magnetic  field  distribution  at  the  boundaries.  This, 
in  turn  can  be  evaluated  from  the  total  currents.  The  following 
relations  are  used: 


7  i  <c  V 


*  “o 


J 

z 


dz 


J 


r 


From  Eq. 7  it  follows  that  J  B 


—  (B2/2h  ) 


(6) 

(7) 
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R  !2<'> 


Thrust 


*»  p 

Jr 

o  u 


o  z^r) 


Bg  2nr  dz  dr 


r  B2i 

Be 

2 

L  °J 


z2(r> 


00 


Bq  can  be  found  by  integrating  (Eq.  6) 

v 


Bg(r,  z)  =  ~  J  s  Jz(s,  z)  ds 


Equation  8  shows  that  thrust  can  be  evaluated  in  terms  of 
magnetic  field  and  Eq.  9  shows  that  magnetic  field  depends  only  upon 
axial  current.  If: 

a)  Current  leaves  cylindrical  anode  of  radius  R 

A 

b)  current  enters  circular  cathode  of  radius.  R  with 

'  c 

uniform  current  density 


^  I  ,  R. 
Thrust  -  (|  +  In  f) 

c 


3.3  Torque  Due  to  (r  J  B  *  r  J  B  ) 

z  r  r  z  * 

In  an  axially  symmetric  volume  with  radial  and  axial  currents 
and  magnetic  fields,  there  will  be  a  torque  which  occurs  when  the  cur¬ 
rent  crosses  the  magnetic  field.  To  evaluate  this  torque,  introduce 
the  vector  potential  Aq. 


f3 


5470-Final 


3-3 


•  9  „  n  , 

Br  ■  -  aT  :  h  mr  a?  (r  V 

The  quantity  (r  Ag)  is  constant  along  a  magnetic  field  line, 
torque  per  unit  volume  is  given  by: 

[  a  (r  A„)  a  (r  Aq)1 

rJB  -  rJB  *  -  J  - - — —  +  J  - r— 5- 

z  r  r  z  r  dr  z  az  ! 

L  -» 


For  the  axially  symmetric  case., 


a  (r  Ag) 


=  0,  hence 


Torque  =  -  J  •  V  (r  Ag) 


-  7  •  (r  As  J) 


where  V  .  J  =  0  has  been  used.  Upon  integration  over  a  volume  R  of 
surface  S,  outward  normal  n, 


(torque)  d(vol) 


-J<- 


r  Afl)  (J  .  n)  dS 
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I 

1 


i 


t 


If  r  Aa  is  constant  at  anode  and  cathode 


Torque  -  I  [(r  Ae)anode  -  (r  A,)^^  ] 


anode  r 


=  I 


U 

cathode 


d  (r  A 


9r 


0)  S  (r  Afl)  *1 

Ldr  +  _^d2j 


(15) 


I  (Magnetic  Flux  Between  Cathode  and  Anode) 

2rr 


Finally  for  a  point  cathode,  and  an  average  axial  field  B  through  a 

z 

circular  anode  of  radius  RA> 

Torque  =  |  B  I  (16) 


3.4  Max  Thrust  Due  to  JJB 

o  r 

In  an  axially  symmetric  volume  where  is  induced  (by  the 
Hall  effect),  the  amount  of  axial  force  cannot  be  larger  than  for  the 
case  of  a  completely  diamagnetic  plasma.  In  this  limiting  case,  the 
Jq  lies  completely  in  the  surface  of  the  volume,  and  no  magnetic 
fields  exist  inside  the  volume.  Compute  the  currents  and  magnetic 
field  as  follows.  Let  B  a  Bq  +  where  BQ  is  the  applied  field  due 
to  external  magnets  and  B^  is  the  field  due  to  induced  currents 
within  the  volume.  Outside  of  the  plasma,  V  x  B^  is  zero,  hence 
Bj.  «  V  Y,  where  Y  is  a  scalar  riuld.  Since  V  .  B^  =  0,  then 

2 

V  Y  =  0  (outside  of  plasma)  (17) 

Since  B  .  n  =  0  at  the  plasma  surface 

=  -  (Bq)  .  n  (at  plasma  surface)  (18) 
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fs~ 


Solve  <Laplace  Eq.  17  with  boundary  condition  Eq.  18.  The  force  is 
given  by 

Force  =  -  J  f  (BQ  +  V  Y)2  n  d  S 

where  n  is  an  outward  normal  and  S  is  the  surface  area  of  the  plasma. 

The  same  result  can  be  obtained,  by  the  following  method. 

Here  we  shall  attempt  to  solve  for  the  surface  Jg  distribution  for  a 
particular  plasma  configuration  shown  in  Fig.  3-1.  We  use  a  0-Ohm's 
law  (where  Yg  is  the  Hall  parameter  for  electrons). 

J0  -  <W> 
and  the  induction  equations 


SB  dB 
V  z 

dz  dr  ^o 


(20-a) 


1  a  S® 

7  dr  <rV  +  bT  *  0  (20-b> 

Y 

In  Eq.  19,  the  coefficient  p~-j  does  not  depend  upon  magnetic 
field,  since  Y  is  proportional  to  J  B  |.  We  assume  J  and  J  are  known, 
based  upon  visual  observations  of  the  operation  of  the  accelerator. 

It  appears  that  =  0,  and  J =  total  current  divided  by  cross  sec¬ 
tional  area  of  the  anode  jet  (2rtR  6).  Thus  Eq.  19  becomes 

a 

JG  ”  '  [ph  Br  <21> 

From  Eq.  21  we  see  that  Br  (and  not  Bz>  is  important  for 
calculating  Jg.  We  can  use  the  integral  solution  of  Eqs.  20-a  and 
20-b  which  assumes  no  magnetic  material  is  present: 
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WIRE  MODEL 


FIG.  3-1  MODELS  FOR  INDUCED  FIELD  ANALYSIS 
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,  t)  ds  dt 


Br  <r*  *>  -  2^11°  (ti  Si  Z‘t)  JS 


where 


(22) 


G(r »  s;  z-t)  s  J* 


r  s  (z-t)  cos0  d 0 


o  j~r2  -  2rs  cos0  +  s2  +  (z-t)2j3/2 

When  applying  Eq.  22,  we  use  Eq.  21  for  the  value  of  Jg  in  the  region 
of  the  anode  jet,  and  the  known  Jg  iu  the  electromagnet.  Thus,  Eq.  21 
becomes : 


r 


r  Ye  *A  “|  j 

J9,Hall  =  L"  Fb|  2  it  R  6j  2*r\ 


J’J  G  (r,  s;  z-t)  Jg^oi]/8**1)*18  dt 


coil 


)  (23) 


'+  JJ  G  (r,s,z-t)JQ^Hall(s,t)  ds  dt 
^anode  jet 


Eq.  23  is  an  integral  equation  for  the  Hall  current  in  the  anode  jet. 

Exact  solutions  are  not  available  for  Eq.  23.  An  approximate 
solution  has  been  obtained  by  lumping  the  distributed  Hall  currents  Jg 
into  concentrated  currents.  We  will  refer  to  this  technique  as  the 
wire  model.  The  Hall  currents  are  replaced  approximately  by  hoop  cur¬ 
rents  in  a  set  of  wires.  The  error  of  the  lumping  has  been  estimated, 
by  changing  the  number  of  wires  per  unit  length  of  the  anode  jet,  and 
shown  to  be  small. 

The  solution  found  was  for  infinite  conductivity,  or  more 
precisely  \  •*  »  where 


\ 


li  Y 
c  e 


B  I  (2ir) 


(24) 


This  is  the  limit  of  a  completely  diamagnetic  plasma. 
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To  solve  the  wire  model,  concentrate  the  current  density 
into  a  set  of  wires  (see  Fig.  3-1).  Thus,  the  wire  for  the  electro¬ 
magnet  coil  has  a  current  I  ,  where  I  is  equal  to  the  number  of  turns 

c  c 

times  the  current  in  one  conductor.  The  wires  which  carry  the  Hall 
currents  in  the  anode  jet  carry  a  current  of  Jg  Hall  •  Az  •  6.  The 
integral  £q.  23  is  thus  replaced  by  a  set  of  simultaneous  algebraic 
equations.  In  the  infinite  X  case,  which  corresponds  to  infinite 
conductivity,  =  0.  Thus, (B^) ^  =  0,  where  j  stands  for  one  of  the 
points  shown  in  the  wire  model  of  Fig.  3-1. 

0  "  G<V  RC-  zj>  Jc  *1  G<V  V  -  *i>  <25> 

i=l 


Some  solutions  of  Eq.  25  are  shown  in  Fig.  3-2.  The  Kernel 
function  G  defined  in  Eq,.  22  can  be  computed  in  terms  of  elliptic 
integrals. 


G(r ,s;t)  -  p  2 


2rst 


i,  ,  it  ,  2H3/2 

j(r+s)  +  t  j 


(2-k2)  E(k)  -  2  C 1— k2)  K(k) 
k2  (  1  -  k2) 


(26) 


2  r  2  2“> 

where  E  and  K  are  elliptic  integrals  and  k  **  4rs/|  (r+s)  +  t  •.  To 

W  aj 

simplify  the  computations,  a  function  G  .  based  upon  a  flat  geon. 

planar 

etry  was  used 


Gplanar  “ 


:  2  2  I  ‘  S 

j  (r-s)  +  tZj  !  (r-i-s) 


(27) 


This  is  a  good  approximation  to  the  Kernel  function  near  the  coil. 

The  electromagnetic  thrust  can  be  computed  either  from  the 
integral  of  the  JfiB  forces  in  the  jet  or  by  the  reaction  on  the  coil. 


Thrust  =  Jg  co^  2itr  dr  dz  = 
coil 


JT  J9,HallBr  2“  dr  dz  <28> 
anode  jet 
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FIG*  3-2  COMPUTED  HALL  WIRE  CURRENTS  AND  THRUSTS 


i  >  i 


' 
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For  the  example  of  Fig.  3-1,  the  integration  was  made  over  the  coil. 
The  thrust  is  approximately  80  percent  of  the  product  of  the  cross- 
section  area  of  the  anode  jet  times  the  magnetic  pressure  which  would 
exist  along  the  centerline  if  there  were  no  Hall  currents.  This  is 
probably  an  upper  limit  of  the  possible  thrust. 


W 


1 

/ 6 / 
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